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I. INTRODUCTION AND SUMMARY 

In this volume there is an introduction to cosmology and CMB physics by Sourdeep, and 
on the inflationary paradigm by Wands. They summarize the synergy between theory and 
observations that has produced spectacular advances in our understanding of the universe 
in the last decades. The emergence of a "concordance model" is a remarkable success 
of cosmology and the theory of General Relativity in which the current paradigm relies. 
However, the widely accepted Hot Big Bang scenario, regarded as the "standard model 
of cosmology", contains important limitations, already manifest in its name. The model 
encompasses a phase in the very early universe in which the density of matter and the space- 
time curvature grow unboundedly, blowing up at the big bang singularity. The big bang is 
not a prediction, but the result of applying the theory beyond its domain of validity. When 
the energy density and curvature approaches the Planck scale, the predictions of General 
Relativity are unreliable; the quantum aspects of the gravitational degrees of freedom are 
expected to dominate in that regime. This chapter provides a possible quantum gravity 
extension of the well established cosmological model from the perspective of loop quantum 
gravity. 

Loop quantum cosmology (LQC) arises from the application of principles of loop quan- 
tum gravity (LQG) [1] to cosmology. The goal is to quantize the sector of General Relativity 
containing the symmetries of cosmological space-times, by following the physical ideas and 
mathematical tools underlying LQG, presented in detail in the chapter by Sahlmann. Re- 
stricting attention to cosmology presents several advantages. The existence of underlying 
symmetries largely simplifies technical issues, and allows to overcome mathematics difficul- 
ties that are hard to handle in more generic situations. Yet, the structure is rich enough 
to contain deep conceptual issues in quantum gravity: What happens with space and time 
when matter density and curvature reach the Planck scale. Does the big bang singularity 
persist? What is the meaning of time in the Planck era? How do classical General Relativity 
and a smooth space-time description arise in the low energy regime? What is the scale at 
which quantum gravity effects become subdominant? Does quantum gravity have anything 
to contribute to the origin of cosmic structures and to the infiationary scenario? On the 
other hand, the astonishing advances in theoretical and observational Cosmology in the last 
years have been able to relate observations with theories of the very early universe. Cosmol- 
ogy then offers an interesting arena in which quantum gravity can make contact with other 
theories such as inflation, and probably provides the most promising avenue to confront 
quantum gravity ideas with observations. 

But the restriction to cosmological settings also leads to important limitations. In princi- 
ple, it is not guaranteed that the result of quantizing a symmetry reduced sector of general 
relativity will reproduce the same physics as the restriction of a full quantum gravity theory 
to symmetric scenarios. Symmetry reduction often entails a drastic simplification, and one 
may loose important features of the theory by restricting the symmetry prior to quantiza- 
tion. However, it has been extremely useful in several areas of physics, when the complexity 
of the problem under consideration made it difficult to find solutions without introducing 
additional inputs. The Oppenheimer-Snyder model of black hole formation, or the Dirac 
quantization of the hydrogen atom are examples that were able to encode the key physical 
ingredients of the problem, in spite of the severe symmetry reduction. Quantum cosmology 
may well be another example, if it is constructed choosing carefully the key ingredients from 
full quantum gravity. It is likely that predictions from quantum cosmology will not agree in 
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every detail with those obtained from full quantum gravity applied to cosmological scenar- 
ios, but we expect it to capture the main aspects of the complete theory. As in the previous 
examples, quantum cosmology can provide valuable information about the correct way to 
quantize gravity, and be as useful as the hydrogen atom has been for quantum mechanics. 

This chapter provides a brief and pedagogical summary of the advances in Loop 
Quantum Cosmology, with some emphasis on recent results. They can be divided into 
three parts, which are in one-to-one correspondence to the three sections in which the 
chapter is divided: 1) Quantization of cosmological space-times; 2) Inhomogeneous pertur- 
bations in LQC; 3) LQC extension of the inflationary scenario. In the remainder of this 
introduction we summarize the content of each of these sections and provide a global picture. 

1) Quantization of cosmological space-times. General Relativity is a totally con- 
strained theory, in the sense that the full Hamiltonian generating dynamics is required to 
vanish. Something similar happens in classical electromagnetism, where part of the Hamilto- 
nian, the piece that generates gauge transformations, is a constraint. In General Relativity 
the constraint turns out to be the full Hamiltonian, reflecting the background independence 
of the theory. Dirac provided the conceptual framework to quantize constrained systems. 
At the quantum level, physical states have to be annihilated by the operator corresponding 
to the classical Hamiltonian, C ^ = 0, and all the physics has to be extracted from this 
equation. The quantum state ^ is the wave function of the physical fields, including the 
gravitational field itself, and classical quantities such as the metric, energy density, curva- 
ture tensor, are represented by quantum operators on the physical Hilbert space Tiphy it 
belongs to. The non-trivial mathematical problem is to make sense and solve the quantum 
constraint equation, and the underlying cosmological symmetries largely facilitate this task. 

The next conceptual issue is to obtain the familiar time evolution that we normally 
use in physics from this time-less or 'frozen' formalism. At the quantum level we do not 
have a classical metric teUing us what are the time-like directions in the manifold, and all 
what we have is a probability-distribution of different metrics. A useful strategy has 
been to follow a 'relational-time' approach, in which one of the physical variables of the 
problem plays the role of time, and the rest evolve with respect to it. By using a massless 
scalar field as this internal time, it is possible to construct the Hilbert space of physical 
states satisfying the quantum constraints, and a precise mathematical framework has been 
developed to to study the resulting quantum geometry [2]. It has been shown that all the 
operators representing physical quantities such as the energy density, space-time curvature, 
etc, remain bounded on the physical Hilbert space, even in the deep Planck regime. This is 
the mathematical sense in which the singularity is resolved in LQC. The physical picture 
that emerges from the abstract formalism is the following. When the energy density of 
the universe is comparable to the Planck energy density, the quantum properties of space- 
time geometry become important and dominate. A sort of quantum repulsive degeneracy 
force appears at such extreme densities, precludes the universe to continue contracting, and 
forces the quantum space-time to expand again once the maximum energy density has been 
attained, replacing the big bang singularity by a quantum bounce. This maximum energy 
density is proportional to h~^, similar to the finite energy of the ground state of the hydrogen 
atom that avoids the collapse of the positron and electron as a consequence of the Heisenberg 
principle. When the energy density and curvature become smaller than approximately one 
percent of the Planck scale, the quantum effects of gravity become rapidly negligible and 
classical General Relativity provides an excellent approximation. The resulting quantum 
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dynamics has been analysed in detail and has provided important insights on the behaviour 
of physics in the Planck regime. The ability of incorporating non-perturbative quantum 
corrections that are able to completely dominate the evolution in the Planck regime and 
dilute the big bang singularity and, at the same time, to disappear in the low energy regime 
to find agreement with the classical description, is a highly non trivial result of LQC. 

Remarkably, some global aspects of the evolution of the quantum geometry can be en- 
coded in simple effective equation. Those equations provide a smooth space-time metric that 
approximates the full quantum evolution of the quantum space-time. They have similar form 
to the equations arising in General Relativity, but include new terms, proportional to h, that 
make the effective trajectory to depart from the classical one around the Planck era. The 
effective dynamics provides an excellent approximation of the quantum evolution, even at 
Planckian densities, provided the quantum state is chosen to be highly peaked in a classical 
trajectory in the low energy regime where General Relativity provides a good approximation. 

2) Inhomogeneous perturbations in quantum cosmology. As emphasized in the 
chapters by Sourdeep and Wands, the theory of inhomogeneous perturbations (of matter and 
gravitational degrees of freedom) propagating in classical cosmological space-times has been 
a key mathematical tool in modern cosmological research. One of the deepest insights in 
cosmology is the idea that the cosmic structures (galaxy clusters, super-clusters, etc) that we 
see today were originated in the very early universe by a process of amplification of quantum 
fluctuation by the cosmological expansion, as explained in the context of cosmic inflation in 
the chapter by Wands. In the inflationary scenario, this occurs when the energy density in 
the universe was close to the GUT scale, (lO^^GcV)^, around 12 order of magnitude below 
the Planck energy density. Quantum gravity effects of the background space-time metric are 
subdominant at those scales, and the theory of quantized fields propagating in a classical 
background appears to be the appropriate mathematical framework to work out physical 
predictions. However, earlier in the evolution of universe, when the curvature and energy 
density are close to the Planck scale, quantum gravity effects are expected to be important, 
and they should not be ignored. To have a complete picture of the evolution of cosmic 
inhomogeneities that encompasses the Planck regime we need to learn how quantum fields 
propagate on a quantum cosmological space-time [3, 4]. The goal of the second section of 
this chapter is to review the construction of such a theory. 

The detailed description of quantum cosmologies provided by LQC is the suitable arena. 
The construction of QFT on quantum cosmologies follows closely the guiding principle be- 
hind LQG: first carry out a truncation of the classical theory adapted to the given physical 
problem, and then quantize by using LQG techniques. The sector of the classical theory of 
interest is extended in this part to cosmological background plus first order inhomogeneous 
perturbations on it. 

The resulting framework originates from first principles, under the assumption that 
inhomogeneous behave as test fields on the quantum geometry, and it should provide a 
bridge between quantum gravity and QFT on curved space-times. Therefore, it is suitable 
to face important conceptual questions such as: What are the concrete approximations 
under which the familiar quantum field theory (QFT) in classical space-times arises from 
this more complete description? What are the precise aspects of the quantum geometry 
that are 'seen' by the quantum fields propagating on it? Does the resulting QFT make sense 
for trans-Planckian modes? These issues will be discussed with some detail in Section 11. In 
Section III, this framework is applied to the study of gauge invariant cosmic perturbations 
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and phenomenological consequences are worked out. 

3) LQC extension of the inflationary scenario. The inflationary scenario occupies 
the leading position in accounting for the origin of the cosmic inhomogeneities observed in 
the Cosmic Microwave Background (CMB) and large scale structure. This success is mainly 
rooted in the economy of assumptions, the elegant mechanism that originates the cosmic 
inhomogeneities from vacuum quantum fluctuations, a subtle interplay between quantum 
mechanics and classical gravitation, and particularly the non-trivial agreement with obser- 
vations. Inflation is however an effective theory, and it is expected that a more fundamental 
theory will complete it. Examples of open questions that the more complete theory should 
answer are: What is the nature of the scalar inflaton field? Is there a single or several 
fields, like in multi-field models? What is the specific shape of the infiaton potential? These 
questions originate in particle physics, and unfortunately at these stages LQC does not have 
much to contribute. There are, in addition, important issues related to gravitation: What 
is the evolution of the space-time before infiation? In General Relativity the big bang sin- 
gularity is unavoidable in inflationary scenarios [5]. Is there a quantum gravity scenario 
in which the singularity is resolved and in which the evolution flnds an inflationary phase 
compatible with observations generically, i.e. without a fine-tuning of its parameters! Such 
a scenario would allow to extend the inflationary space-times all the way back to the Planck 
era. Moreover, one could then use the quantum theory of cosmological perturbation on quan- 
tum space-times described in section III, to extend the analysis of cosmic inhomogeneities 
to include Planck scale physics. 

Section IV will review the arguments showing that such an extension is possible in 
LQC, where one can construct a conceptual completion of the inflationary theory from the 
quantum gravity point of view, in which Planck scale physics can be included in the study 
of cosmological perturbations. The importance of this extension goes, however, beyond the 
conceptual domain and may open a window for phenomenological consequences. 

To summarize, this chapter will review recent advances in the completion of the quantiza- 
tion program underlying LQG when restricted to the cosmological sector. We shall explore 
how the singularity of the homogeneous background is avoided, and how the abstract theo- 
retical framework can descend down to make contact with phenomenology. Although many 
open issues still remain, at the present time there is a solid body of knowledge, based on a 
rigorous mathematical framework. These combine with analytical and numerical techniques, 
and provide an avenue from the big bang singularity resolution to concrete observation of 
the CMB and galaxy distributions. 

Due to space restrictions, there are some topics that we shall not cover in this chapter, 
such as the path integral formulation and its relation with spin foams [6], spin foam cos- 
mology [7], the Gowdy models [8-12], nor numerical issues [13]. We do not provide either 
a review of all the existing ideas to study LQC effects on cosmic perturbations. See [14-23] 
for different approaches to that problem. Further information can be found in the reviews 
[24], [25], [26] and [27]. 

Our convention for the metric signature is — we set c = 1 but keep G and h explicit 
in our expressions, to emphasize gravitational and quantum effects. When numerical values 
are shown, we use Planck units. 
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II. QUANTIZATION OF COSMOLOGICAL BACKGROUNDS 

In this section we shall consider the quantum theory of the homogeneous background 
within the context of Loop Quantum Cosmology. First we shall discuss what it means 
for a cosmological model to be quantized, or to use the standard nomenclature, to define 
a quantum cosmology. Just as with the quantization of any mechanical system such as 
the hydrogen atom, the first step is to cast the model to be quantized in a Hamiltonian 
language. That is, one has to identify configuration variables and their corresponding 
momenta pj, with the property that the Poisson bracket is {q\Pj} — Sj. The next step 
in the quantization process is to find a Hilbert space H and operators and pj satisfying 
[q^,Pj] = ihSj. Then one has to define an operator H corresponding to the Hamiltonian 
(and to other physically relevant observables) , in order to define dynamics through the 
Schrodinger equation: —ihdt^ = H'^. 

In the case where the classical system under consideration is a totally constrained sys- 
tem, instead of a Hamiltonian H defining dynamics, both the classical description and the 
corresponding quantization are more subtle. Here the dynamical variables are subject to a 
constraint C{q,p) = 0. Furthermore, there is no Hamiltonian defining dynamics, and the 
canonical transformations generated by the constraint C are interpreted as gauge. That 
is, points on the phase space connected by a canonical transformations generated by the 
constraint are physically equivalent. Thus, the curve on phase space made out of all the 
physically equivalent points represents a gauge orbit and can be identified with a point on 
the true, physical phase space. Observables will be those functions f{q,p) that are constant 
along the gauge orbits (i.e. satisfying {/, C} = 0). Since there is no true dynamics, the 
system is said to posses a frozen dynamics. A natural question is whether one can extract 
some 'dynamics' from the frozen formalism. In some cases, one can use one of the variables 
(or an appropriately selected function) as an internal time T{q,p), with respect to which the 
gauge orbit can be described in terms of a relational dynamics (that is, where the 'dynamics' 
is described by correlations between the variable T and the rest of the variables) . 

Let us now review the quantization process when we have a totally constrained system. 
The first step is to define a kinematical Hilbert space "Hkin- This space serves as an arena for 
the implementation of the constraint, that is now required to be represented as a self-adjoint 
operator C on "Hkin- Not all states in the kinematical Hilbert space are regarded as physical. 
The condition that selects those physical states was put forward by Dirac and has the form, 

C-^phy = 0. (2.1) 

Once one has found the physical states ^^phy (that might belong to "Hkin or not), one needs 
to specify an inner product (-l^hy in order to construct "Hphy, the physical Hilbert space. 
Physical observables will be operators F that leave the space of physical states invariant. 
This translates into the condition [F, C] = 0. In some cases, when there is an internal time 
variable T, one can recast the Dirac condition (2.1) as an 'evolution' equation where T plays 
the role of time, as in the Schrodinger equation. 

One interesting feature of the simplest cosmological models is that they are totally con- 
strained systems, so the general framework we have outlined is applicable. Even more, one 
can complete the quantization program and obtain a complete physical description where a 
massless scalar field plays the role of internal relational time. One can then pose physical 
questions pertaining to observables of cosmological interest, such as the Hubble parameter 
and curvature scalars. Interestingly, for the simplest models, one can indeed find two dif- 
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ferent, inequivalent, quantizations. The first one corresponds to the so-called Wheeler-De 
Witt (WDW) quantization that was put forward by De Witt and Misner in the 60 's. The 
second quantization corresponds precisely to the one we shall here consider in detail, known 
as loop quantum cosmology. As we shall describe in more detail later, the basic difference 
between these two programs corresponds to the choice of kinematical Hilbert space "Hkin- 
The choice made by De Witt and others was, in a sense, the most natural one, resembling the 
Schrodinger quantum mechanics that has been so useful to describe many physical systems. 
On the other hand, the choice one makes in LQC is somewhat exotic from the perspective 
of standard quantum mechanics, but is selected when the underlying symmetries pertinent 
to the gravitational field are seriously taken into account. 

The second and physically most important difference between these two representations 
is that their predictions regarding the fate of the classical singularity are radically different. 
While the WDW theory predicts that the singularity remains, as defined by the behavior of 
the expectation values of physically relevant operators such as energy density, in the case of 
LQC the singularity is generically avoided. Instead of a big bang (or big crunch) one has 
a bounce connecting a contracting branch with an expanding one; the energy density and 
curvature scalars are bounded from above, so that physics is well defined throughout the 
intrinsic dynamical evolution of the quantum state describing the universe. 

Let us now briefly describe the structure of the remainder of this section. In the first 
part, we study in detail the k=(} FLRW model with vanishing cosmological constant, and 
discuss some of its main features. In the second part we discuss other models. The first one 
we consider is the closed k—1 model also without a cosmological constant. Next, we briefly 
discuss k=0 FLRW models with a cosmological constant and some anisotropic models. In 
the third part, we introduce the so called effective equations. We give a brief introduction 
to the subject and discuss in detail the case of k=0 FLRW model. Next we consider the 
k—1 case, followed by a discussion of anisotropic effective space-times, including the Bianchi 
I, II and IX models. 



A. A; = FLRW, singularity resolution 

The simplest model that one can consider is a k=0 homogeneous and isotropic FLRW 
cosmological model foliated by 3-manifolds S that are topologically M^. In order to find a 
Hamiltonian description for the model, we have to start with an action principle. Due to 
homogeneity, the action is not well defined unless one introduces and fixes a fiducial cell V. 
This will play the role of a co-moving volume. We can introduce a fiat fiducial metric Qab 
on with respect to which the coordinate volume of V is K = ^/qd'^x. Without loss of 

generality, in what follows we shall set = 1. The fiat FLRW spacetime is described by 
the metric 

ds' = -NMt'' + a(i)Mx2 (2.2) 

where N is the lapse function, q dx^ is the ffat fiducial metric, and a is the scale factor 
of the universe. Now, the action principle is 

J 3, 
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with R the scalar curvature of the spacctimc. The gravitational part of the phase space 
consists of a and its conjugate momenta that is found to be: 

In this simplest model, the matter we shall consider is a homogeneous massless scalar field 
0. The action for such a field is: 



'matt 



di 



N 



From this, the momenta p(^) associated to the scalar field is p(^) = ^j^, and the Hamiltonian 
constraint that defines the 'dynamics' is then. 

To summarize, the phase space is four dimensional with coordinates (a. Pa] (l),P{(f,)), satisfying 
{a, Pa} = 1 and {(f>,P(4,)} = 1. In the standard Wheeler-De Witt approach, the next step 
is to consider the kinematical Hilbert space to consist of 'wavefunctions' ^wdw = 
of the 'configuration' variables (a, 0). In this case, the operators are represented in the 
usual fashion, as: a ■ \l'(a, </)) = a\E'(a, 0) and Pa = —iHda'^{a, 0), and similarly for the other 
variables. Then, one promotes the constraint (2.3) to an operator, and finds solutions to the 
Dirac condition (2.1). This has been described in detail in [28, 29]. 

In order to define the corresponding phase space in loop quantum cosmology, we need to 
follow some more steps. The first one is that one needs to introduce a new set of variables 
for the gravitational degrees of freedom. As explained in Sahlmann's contribution to this 
volume, loop quantum gravity, and consequently LQC is based in a connection A and its 
corresponding momenta E, a generalization of the magnetic potential and electric field of 
electromagnetism. Let us then write the phase space in terms of these so-called Ashtekar- 
Barbero variables. First, introduce a fiducial triad e" and co-triad uj^ compatible with Qab- 
The conjugate phase space variables are the SU (2) connection — + ^K^ and the 
densitized triad satisfying 

{Al{x),E'^{y)} = 87rG^S'j]S'{x,y). (2.4) 

Here is the spin connection measuring the intrinsic curvature (which vanishes in the A; = 
model), 7 is the Barbero-Immirzi parameter and i^* is the extrinsic curvature 1-form related 
to the extrinsic curvature Kab as — e^Kab, with e" the un-densitized triad. Due to the 
underlying symmetries of the homogeneous isotropic spacetimes we are considering, these 
variables can be written as [2] 

Al = cul ; El =p^qe1. (2.5) 

Thus, the dynamical variables in the isotropic cosmological regime are p and c. The rela- 
tionship between the 'triad' p and the scale factor is, 

bl = ■ (2.6) 
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The connection component gets related to the rate of change of scale factor as 

holding only for the physical solutions of General Relativity (GR). The gravitational part 
of the phase space is characterized by the conjugate variables c and p satisfying: 

{CP}-"^. (2.8) 

and the complete phase space has coordinates (c,p; The dynamics thus fund in 

the Hamiltonian language is completely equivalent to the standard description based in 
Einstein's equations. To see that, one can find the Hubble parameter H = a/a = p/{2p) 
by computing p = {p,C}, where C is now written in terms of the variables (c,p; (f),p((f,) ) (see 
(2.12) below). Prom there one can write the standard Friedman equation: H"^ — p, with 

Let us now consider the issue of quantization. As previously discussed, the choice of 
kinematical Hilbert space in LQC is different from the WDW case. That is, we do not expect 
to represent c and p as multiplication and derivation, for example. The idea instead is to 
construct a quantum theory that is closest to the quantization used in loop quantum gravity, 
as discussed for instance in [1]. This means in particular a different choice of kinematical 
Hilbert space. Recently this polymeric quantization for cosmological models has been shown 
to be unique when invariance under diffeomorphisms is imposed [31] (in complete analogy 
with the corresponding results in full LQG [32, 33]). The new strategy is the following. 
Instead of re-writing the Hamiltonian constraint (2.3) in terms of the (c,p) variables, one 
starts with the full expression of the Hamiltonian constraint, in terms of variables A and 
E. Then, one uses the simplification given by Eq.(2.5). As it turns out, the choice of the 
polymeric Hilbert space as the kinematical arena for the implementation of the constraint 
-following the LQG route to quantization- has the important feature that it does not admit 
the c operator. That is, only exponential functions of the gravitational connection c such as 

h[^^^ = cos(Ac c/2)I + 2 sin(Ac c/2)rfc (2.9) 

become well defined. These objects have the geometrical interpretation of being the 
'holonomies', or parallel transports of the connection A. These functions generate an al- 
gebra of so-called almost periodic functions whose elements are of the form exp(iAcc/2). 

The resulting kinematical Hilbert space is then L^(MBohr, d/iBohr), a space of square in- 
tegrable functions on the Bohr compactification of the real line. Beside the name of the 
space, it is straightforward to understand the nature of this space. For instance, the eigen- 
states of p, labelled by |//), satisfy (/Xi|/X2) = S^ll,^i2■ This is to be contrasted with the usual 
Schrodinger representation where, instead of the Kronecker delta one has the Dirac delta. 

In particular, this eigenstates are normalized and constitute a basis for the kinematical 
Hilbert space H^oiy This constitutes the main difference from the standard Schrodinger 
representation where the eigenstates of momentum p = /x are not normalized and 
satisfy = 5{ijl,v). Note also that this plane waves states are not a basis for the 

L^(M, dx) Hilbert space. 

There exists an important result in mathematical physics stating that for a finite di- 
mensional phase space, the Schrodinger Hilbert space is the only choice of representation 
of the canonical commutation relations, satisfying some regularity conditions. This result 
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goes under the name of the Stone- Von Neumann uniqueness theorem [30]. Thus, one could 
have imagined that, since the system has a finite number of degrees of freedom, both the 
WDW and the LQC representations should be equivalent. However, that expectation is 
not realized. The polymeric representation used in LQC and the standard one are unitarily 
inequivalent. This is due to a crucial property of the LQC operators, implying that the 
polymer quantum mechanics does not posses some of the regularity conditions that go into 
the hypothesis of the Stone- Von Neumann theorem. To explore those propertie further, let 
us consider the action of the two fundamental operators on the eigenstates 



= — ; exp(iAec/2)|/x) = |/x + Ae) (2.10) 

Note that the 'displacement' operator exp(iAcc/2) is not continuous when Ac — > 0, since the 
states and |yU + Ac) are always orthogonal to each other, for all values of Ac > 0. Also note 
that a basis of the polymer Hilbert space is uncountable as the label n for the eigenstates 
can take any value in the real line. 

In order to obtain the quantum constraint the key step is to rewrite the classical gravi- 
tational constraint with field strength F^j, as. 



f E"' E 

Cg.av = -7"' ei.fc Kb (2.11) 



Further, one writes the field strength F*^ in terms of holonomies and triads and then quantize 
(where we have chosen = 1). The matter part of the constraint is quantized in the regular 
way, where the Schrodinger representation is used. A further simplification is to choose 
N — — V from the very beginning. If we rewrite the line element with this choice we 
have ds^ = — a^dr -|- a^dx^, for which the classical constraint reads. 

With this choice, the gravitational constraint has then the form, 

Cgrav = -7-V\e«e$F„^, (2.13) 

The field strength can be classically written in terms of a trace of holonomies over a square 
loop Dy, considered over a face of the elementary cell, with its area shrinking to zero: 



with 

^S;.^ = hf''^hf^\h^^Y\h';^''Y' ■ (2.15) 

Since the underlying geometry in the quantum theory resulting from LQC is discrete, the 
loop Oij can be shrunk at most to the area which is given by the minimum eigenvalue of 
the area operator in LQG: A = kip^ with k of order one. Note that it has been standard in 
the LQC literature to choose k = 2\/37r7 [2], but it can also be taken as a parameter to be 
determined [29]. The area of the loop with respect to the physical metric is A^|p|. Requiring 
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the classical area of the loop to have the quantiim area gap as given by LQG, we are led 
to set Ac = a/A/|p|. Since Ac is now a function of triad, the action of exp(iAc(p)c) becomes 
complicated on the states in triad {/j,) basis. However, its action in volume {u) basis is very 
simple: it drags the state by a unit affinc parameter. 

It is then convenient to introduce the variable b := j^jt^; such that AcC = Abb where 

Ab := "n/A is the new affine parameter. Note that b is conjugate variable to u, satisfying 
h{h, u} — 2, where u labels the eigenstates of the volume operator 

V\u) = 27r4i7kl W) ■ (2.16) 
The action of the exponential operator then becomes very simple: 

exp{iXc/2)\u) = exp{i\^/2)\u) = |z/ + Ab). (2.17) 

In what follows we shall only consider Ab a constant and shall only denote it by A. Further, all 
of the identities used to write classical constraint in terms of holonomies remain unaffected 
and the quantum constraint operator on wave functions \&(z/, 0) oi u and is obtained 

.^9 ? / ,v ^ ^ I I sin Ab , , sin Ab ~ , 

dl ^{u, 4>) = SttG \v\ -j— \iy\ *(i/, <j>) (2.18) 

Writing out the explicit action of operators sin Ab, (2.18) simplifies to: 

o2 ? / i\ o ^ sin Ab sin Ab ~ , ^ , 
aJ^'(j/,0) = 3TTGiy^—iy^—-^{iy,(l)) 



3tiG 

V 



+ 2A)§(i/ + 4A) - 2v^{v) + (i/ - 2A)^(i/ - 4A) ] 



4A2 

=: e(,)^(z/,0). (2.19) 

The geometrical part, ©(i/), of the constraint is a difference operator in steps of 4A. 

C^{v)^{v + 4A) + CH{v) + C-*(i/ - 4A) = C^att "^{v) (2.20) 

where and C° are functions of [34]. Note that the equivalent of the Wheeler-De Witt 
equation is now a difference equation in the geometrical variable, instead of a differential 
equation. 

Then, physical states correspond to solutions to the quantum constraint (2.19), but they 
should also belong to the positive frequency part of the Hamiltonian constraint, and satisfy 
the 'Schrodinger equation', 

-in9^*(i/,0) = V©*(i/,0). (2.21) 

Furthermore, they should be symmetric under v — > —v and have finite norms under the 
inner product, 

(*l|*2) = *l(^>o) |i^r'*2(i^,0o) ■ (2.22) 

V 

where the constant (^o is arbitrary. As discussed above, these physical states can be inter- 
preted as being solutions to 'evolution equations' with respect to the internal time 0. One 
can indeed define 
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The next step is to define relational observables that will have a clear interpretation in 
terms of 0. For instance, one can define the operator V^^, as the operator corresponding to 
the volume V when the scalar field takes the value 4>q. One can indeed define such Heisenberg 
operators by the standard prescription: 

VU, ■ *phy(i/, 0) := e^^/^^^-^o) v^g-.^/e(0-<Ao) ^^^^^^^ (2.23) 

where V is the standard Schrodinger operator (acting by multiplication in this case). In this 
manner one can define operators corresponding to matter energy density p<^Q and curvature 
scalars, all with a clear interpretation as being defined at 'time 0o'. 

As it turns out, one can perform a Fourier transform into the conjugate variable to u, and 
the resulting quantum constraint, a differential equation, can be solved exactly [29]. This 
allows one to have closed expressions for the expectation values of the Heisenberg operators. 
Let us now describe this solvable model within LQC. 

Solvable loop quantum cosmology (SLQC). We now wish to work in the b representation be- 
cause the geometrical part of the quantum constraint will also become a differential operator. 
Since \^(z/, 0) have support on the 'lattice' v — 4nA, and since b is canonically conjugate to 
i/, their Fourier transforms ^(b, 0) have support on the continuous interval (0, tt/A): 



*(b,0):= y e^'^^(z/,0); so that ^(z/,0) = - / dbe-^^*^ ^(b, 0) . (2.24) 



From the form (2.18) of the constraint it is obvious that it would be a second order differential 
operator in the b-representation. Let us set x(^)0) = Then, on x(b,0), the 

constraint (2.18) becomes 

dl x(b, 0) = ' x(b, 0) (2.25) 

with a — \JV1'kG. Note however, that we did not arrive at (2.25) simply by replacing b in 
the expression of the classical constraint by sinAb/A as is often done. Rather, (2.25) results 
directly from the 'improved' LQC constraint if one begins with a harmonic time coordinate 
already in the classical theory. 

To simplify the constraint further, let us set 

X = - ln(tan ^), or b = ^ tan"^ (e^^ (2.26) 
a 2 A 

so X ranges (— oo, oo). Then (2.19) becomes just the Klein-Gordon equation 

dl x{x, 0) = dl x{x, 0) =: -e x{x, 0) . (2.27) 
The physical Hilbert space is given by positive frequency solutions to (2.27), i.e. satisfy 

-idMx,<f>)^Vex{x,<l>). (2.28) 



We can again express the solutions in terms of their initial data and decompose them into 
left and right moving modes x{x,4>) = Xl{x+) + Xr{x-)- The physical states that we 
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shall consider are positive frequency solutions of (2.27). Since there are no fermions in the 
model, the orientations of the triad are indistinguishable and 0) satisfy the symmetry 
requirement x{~^^ 0) = ~x(^) 0)- Thus, we can write x(a;, 0) = (F(a;+) — F(a;_))/\/2, where 
F is an arbitrary 'positive frequency solution'. To be precise, F{x) is a positive momentum 
function, i.e. with a Fourier transform that has support on the positive axis. With such a 
choice, the solution to the constraint equation become of positive frequency. The physical 
inner product is given by. 



{Xl,X2, 



phy 



-t I [xi{x,(l))d^X2{x,(j)) - {d^xi{x,(j)))x2{x,(j))]dx (2.29) 

/oo 

[^,F^{x+)F2{x+) - d,F,{x_)F,{x_)] dx . (2.30) 
-oo 



The action of the operator on physical states is then: P(<p)X — ~^^5(^X = \f~^ X- 
We can now compute the expectation values and fluctuations of fundamental operator such 
as V^l^o, and p(<^). For any state on the physical Hilbert space the expectation value of the 
volume operator at 'time 0' is given by 



— (X, "^'Ux)phy = 27r74i(x, |i^lx)phy 
where l^^l is the absolute value operator obtained from 

2A 



V — 



a 



cosh.(ax)idx ■ 



(2.31) 



(2.32) 



Using the inner product (2.30) the expectation value of |i>| is given by 



(x, l^^lx) 



phy 



2A 



a 



[d,F{x+){uF{x+)) - d,F{x_){-uF{x_))] dx 

[dxF{x+) cosh.{ax)dxF{x+) + d^F^x^) cosh.{ax)dxF {x_)] dx 



4A '"^ 
a 



dF 



dx 



cosh(Q;(a; — (p)) dx . 



(2.33) 



Prom these expressions one can find the expectation value of certain relational (Heisenberg) 
operators. For instance, the expectation value for the volume operator, at time 0, takes the 
form, 

(\/)^ = \4e"*^ + V^e-""^ , (2.34) 
with, V± constants that depend on the details of the initial (normalized) wave- function: 



± 



47r74iA 



dF 



dx 



dx 



(2.35) 



^From (2.34), it follows that the expectation value of the volume V^l^ is large at both very 
early and late times and has a non-zero global minimum 



K.in = 2(i4y_)V2. 
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The bounce occurs at time 

Around (j) = cf)^ , the expectation value of the volume {V)tp is symmetric. Thus we see that 
all states undergo a big bounce that replaces the big bang (in which the volume goes to zero 
as — )■ ±00). Note: In the case of the WDW quantization, the expected volume reaches 
zero as — >■ ±00, so in this sense one still reaches the singularity. 

Another important observable to consider is the energy density p|<^o. Interestingly, this 
quantity possesses an absolute upper bound on the physical Hilbert space. Let us now see 
how this bound for energy density arises. Fix any state x(a;, 0) = (l/\/2)(F(a;+) — F{x-)) 
in T-Lphy Let us work in the Schrodinger picture at a fixed instant of time, say (pQ. Then, it 
follows that p— {p\ct>o) is given by [29], 



r dx\d,F[ ' 



/^dx|9a;Fp cosh(Q;x) 



2 



(2.36) 



where the integrals are performed at = 0o- Since cosh(aa;) > 1, it follows that the ratio 
of the the two integrals is bounded above by 1. This immediately implies that there is an 
absolute bound given by 

3 1 

{P)<P < Pmax with pmax := ^^^2Q ^2 ' (^■^'^) 

It is interesting to note that this quantity depends inversely with the loop quantum geometry 
scale A. Thus, in the limit A — > 0, where we expect to recover the WDW theory, the density 
becomes unbounded. That is precisely what is found in the complete quantization of the 

WDW theory [29]. 

Using the standard choice for A in LQC, namely A^ = An^/S^ipi, we obtain: 
Pmax = {\/S)/{32TT'^'y^G'^h) ~ 0.41ppi, where we have used the standard choice for 7 coming 
from the black hole entropy computation in LQG where 7 = 0.237. 

In a similar fashion, it is straightforward to see that one can also bound the corresponding 
operator for the 'Hubble parameter operator' Hl^ for physical states. In this case the bound 
takes the form, (H)^ < V (2A7). Note that, just as in the case of energy density, the bound 
on the Hubble parameter is inversely proportional to the loop quantum cosmology scale A. 
In the limit A — > 0, the corresponding quantity becomes unbounded. 

Let us now summarize the main features of the complete quantization of this simple 
cosmological model 

1. The bounce is not restricted to semi-classical states but occurs for states in a dense 
sub-space of the physical Hilbert space. 

2. There exists a supremum of the expectation value for the energy density. This maxi- 
mum allowed density is pmax G'^h). We note that existence of an abso- 
lute maximum of the energy density in this cosmological model implies a non-singular 
evolution, in terms of physical quantities. The singularity is therefore, resolved. 

3. When curvatures become much smaller than the Planck curvature (or for p <^ Pmax) 
the expectation values of the Dirac observables agree with the values obtained from 
classical GR. 
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4. For states which are semi-classical at late times, i.e. those which lead to a large 
classical universe, the backward evolution leads to a quantum bounce in which the 
energy density of the field becomes arbitrarily close to Pmax ~ 0.41ppi. 

5. States that evolve to be scmiclassical at late times, as determined by the dispersion in 
canonically conjugate observables, have to evolve from states that also had scmiclassi- 
cal properties before the bounce (even when there might be asymmetry in their relative 
fluctuations without affecting semiclassicality) [35-37]. Semiclassicality is preserved 
to an amazing degree across the bounce. 

This concludes our discussion of the quantization of the homogeneous background in the 
case that the matter content is a massless scalar field. This is the simplest isotropic model 
and is completely solvable. The question now is how to generalize these results for other 
isotropic and anisotropic models. That will be subject of the next subsection. 

B. Other cosmologies 
1. k=l FLRW 

There are several generalization one might consider away from the k=0, A=0, FLRW 
cosmology. The simplest case is to consider the k—1 FLRW cosmological model [38-40]. 
Even when it is not phenomenologically favored, it is important since it represents a spatially 
closed model that in the classical theory has both an expanding and a contracting phase 
continuously joined by a 'recollapse' point where H = ^ = 0. Therefore, it is an important 
test if one can recover the classical recollapse from the quantum theory. 

The spacetimes under consideration are of the form M = E x R, where E is a topological 
three-sphere S^. It is standard to endow E with a fiducial basis of one- forms 'L'* and vectors 
"e";. The fiducial metric on E is then "g^fc := kij, with kij the Killing- Cartan metric 

on su(2). Here, the fiducial metric "qab is the metric of a three sphere of radius Qq. The 
volume of E with respect to °qab will be denoted by Vq — 27r^ Oq. We also define the quantity 
£o :— Vq^^. It can be written as £q —: i^ao, where the quantity :— {27:'^y^^ will appear in 
many expressions. 

The isotropic and homogeneous connections and triads can be written in terms of the 
fiducial quantities as follows, 

K-yK ; Et^^^V%X- (2.38) 

Here, c is dimension-less and p has dimensions of length. The metric and extrinsic curvature 
can be recovered from the pair (c,p) as foUows, Qab = ^-w°Qab, and "^/Kab = (c — ^) ^-^°Qab- 
Note that the total volume V of the hypersurface E is given by = ^From here, one 

can calculate the curvature of the connection on E as, F^j, = eij^ lo^^lol. The 

only relevant constraint is the Hamiltonian constraint that has the form, 

^^^^^ ^ "8^^ ^ " + ^^-^^^ 

It is convenient to also use the variables [29]: b := c/\p\^/'^ and V = p^/^. The quantity V 
is just the volume of E and b is its canonically conjugate, {b, V} — AnG^. We can then 
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compute the evolution equations of V and b in order to find interesting geometrical scalars. 
Then, 

V = {V, Cgrav} = ^ (bV^ - ^V^^') (2.40) 

from which we can find the standard Friedman equation using the constraint equation C — 

Cgrav + Cmatt ~ and C^M = Vp, WB have := j = ^ P - 

The basic strategy of loop quantization, just as in the A;=0 case, is that the effects of 
quantum geometry are manifested by means of holonomies around closed loops to carry 
information about field strength of the connection. In order to define the quantum theory, 
taking again N = a^, one can work in the u representation and define operators associated 
to curvature and spin connection to arrive to a difference operator Q(k=i) of the form. 



= -e*(i/, 0) + — z/ 



sin^ 



^(z/,(/.), (2.41) 



with K — 27r^i-p\. Numerical solutions of this equation were studied in detail in [38] for 
sharply peaked states, and were shown to posses not only a bounce very close to the critical 
density Pmax, but also a recollapse at a density and volume very close to the classical value. 
Thus, this model provides a very striking example of a quantum gravitational system that 
possesses satisfactory UV and IR behavior. The relative dispersion of V"],^ does increase but 
the increase is very small: For a universe that undergoes a classical recollapse at ^ 1 Mpc, 
a state that nearly saturates the uncertainty bound initially, with uncertainties in and 
V\(f, spread equally, the relative dispersion in V\(f, is still ~ 10~® after some 10^° cycles [38]. 
The expectation values of volume has a quantum bounce which occurs at p = pmax up to the 
correction terms of the order of ^pi/Ko^unce- universes that grow to macroscopic sizes, 
the correction is totally negligible. For example, for a universe which grows to a maximum 
volume of IGpc^, the volume at the bounce is approximately 10^^^£|j. On the other hand, the 
numerical simulations show that one indeed recovers the recollapse with very large precision 
for semiclassical states that reach large volumes [38]. An important lesson that this model 
teaches us is that energy density and curvature are the relevant quantities to define what 
the Planck scale is, and not the size of the universe at the bounce (that, as we have seen, 
can be very large in Planck units). One should also note that, while semiclassical states 
alternate between the Planck scale (UV) and the low density, large volume GR regime (IR), 
states that are 'truly quantum' -or far from semiclassical- might have a bounce at a density 
much lower than Pmax, and not grow to large volumes before recollapse. 

There exists another quantization in which the curvature is not obtained by means of 
closed holonomies, but rather by approximating the connection by open holonomies, as is 
done in anisotropic models with non-trivial curvature [40]. The structure of the constraints 
is different but its quantum solutions have not been explored numerically yet. 

Let us comment on the quantization of the k=-l case. Some early attempts to find such a 
quantization were put forward in [41, 42], but those efforts still suffer from some drawbacks, 
such as absence of essential self-adjointness. A quantization based in open holonomies as in 
[40] is still to be constructed. 
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2. FLRW with A^Q 

The results found for a zero cosmological constant can be generalized to the case of a 
non-zero cosmological constant. For a mass-less scalar field and both signs of the constant, 
we have singularity resolution, in the sense that the big bag/crunch is replaced by a bounce, 
just as in the A = case. For simplicity we shall consider the A < 0, k=0 case, but the 
results can be generalized to k=l as well. The Hamiltonian constraint, for = 1, takes the 
form. 

One can solve the equations of motion and express the dynamics in terms of the scalar field 
as, 

V3|A] cosh[a(0 - 0„)] ^ ■ ^ 

With this, there is a big bang singularity in the past — >■ — oo and a big crunch in the future, 
when — >■ oo. There is a point of recoUapse, when the volume reaches its maximum value 
Knax = (ap(<^))/ (^3|A|), at = 0o, with some resemblance to the k—1 case.The quantum 
constraint takes now the form, 

dl^(u, 0) = -e *(^, 0) - ^^^^ u'^iu, 0) , (2.44) 

with 9 the operator corresponding to the k=0, A=0 case. The operator can be consistently 
defined, and numerically solved [43] to give a picture very similar to the k—1 case with 
vanishing cosmological constant. The big bang/crunch is replaced by a bounce, in such 
a way that a sharply peaked state goes through a series of bounces and recoUapses in an 
almost periodic fashion. 

Let us now consider the A > case. The solution to the classical equations is shghtly 
different from the negative case and takes the form [44, 45] , 

This is qualitatively very different from the previous case. Now, an expanding solution with 
a big bang singularity at the past, — )■ — oo, reaches an infinite volume for a finite value of 0, 
namely when = 0o- Similarly, there are contracting solutions that 'start', for = 0o, with 
an infinite volume and end in a big crunch singularity when oo. At the point = 0o, 
the proper time diverges and the matter density vanishes. One can see that one can actually 
continue the classical evolution past this 'singular' point [45]. In the quantum theory, this 
new behavior manifests itself in the fact that the operator ©a fails to be essentially self- 
adjoint, and one has the freedom of choosing different self-adjoint extensions. Interestingly 
enough, for all of them, the evolution of semiclassical states is almost indistinguishable. 
Evolution is well defined past the point = 0o and the universe recoUapses. As in all 
previous cases, the big bang/crunch singularity is replaced by a bounce. 



3. Anisotropic Cosmologies. 



Isotropic loop quantum cosmology, as we have seen, enjoys a very robust formulation; 
one has complete mathematical control over the quantum theory, one can make physical 



19 



predictions using analytical or numerical tools and can therefore draw conclusions about the 
behavior of a background isotropic quantum geometry. The same is not true for anisotropic 
solutions. While the quantum constraints have been formulated in several cases, one does 
not have full mathematical control regarding their time evolution, and one has not been able 
to solve, even numerically, their dynamical evolutions. In this part we shall summarize the 
formulation of the quantum models as we currently understand them. 

Let us consider the spacetime of the form M = E x M where E is a spatial 3-manifold 
which can be identified by the symmetry group of the chosen model and is endowed with 
a fiducial metric °qab and associated fixed fiducial basis of 1-forms "a;* and vectors °ef . If 
E is non-compact then wc fix a fiducial cell, V, adapted to the fiducial triads with finite 
fiducial volume. We also define Lj which is the length of the iih. side of the cell along °ej 
and V = L1L2L3. We choose for compact E, Lj = V^^^ with i = 1,2, 3. 

Since all of the models in which we are interested are homogeneous and, if we restrict 
ourselves to diagonal metrics, one can fix the gauge in such a way that has 3 independent 
components, c\ and Ef has 3 independent components, pi, 

K=T< ^r = ^\/^X (2.46) 

where pi, in terms of the scale factors a^, are given by \pi\ — LiLjUjak {i ^ j ^ k). Using 
{c\pi) for anisotropic models, the Poisson brackets can be expressed as {c\pj} — 8nGj5j. 
With this choice of variables and gauge fixing, the Gauss and diffeomorphism constraints 
are satisfied and the only constraint is the Hamiltonian constraint 



Ch = I N 



V 



^ ^ Q /I I \ tiu \ ■ I / uu / '^matter 

167rG7^ 



d^x , (2.47) 



with N the lapse function, T^matter = pV and VLab the curvature of the spin connection 
compatible with the triads. 

Using a strategy similar to the isotropic case, the field strength F^^ is given by 

The strategy to choose the corresponding loops is slightly different from the isotropic case. 
We take ii[ = jliLi where fii is a dimensionless parameter and, by previous considerations, 
is equal to jli = Aa/H/ \/\PjPk\ {i j k) [47, 48]. 

For Bianchi II and IX models, this strategy fails because the resulting operator is not 
almost periodic. Therefore, we express the connection in terms of holonomies and then 
use the standard definition of curvature F^^. The operators corresponding to the connection 
are given by [49] 

= ^ . (2.49) 

Note that using this quantization method for fiat FLRW [28] and Bianchi I [47] models, one 
has the same result as the direct quantization of curvature F^j^ (with proper identification of 
the parameters), but for a closed FLRW it leads to a different quantum theory which is more 
compatible with the isotropic limit of Bianchi IX [40, 46, 50]. We call the first method of 
quantization curvature based quantization and the second one connection based quantization. 
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In Bianchi II and Bianchi IX models the terms related to the curvatures, and fi^^,, 
contain some negative powers of Pi which are not well defined operators. To solve this 
problem we use the same idea as Thiemann's strategy, 

|„.|(^-l)/2_ V\Pi\Li , (/i,). I |f/2|^, (2^()) 

where fli is the length of a curve, i E (0, 1) and j e |N is for the representation. Therefore, 
for these three different operators we have three different curve lengths {fi,fi',fl) where fj, 
and jj can be some arbitrary functions of Pi, so for simplicity we can choose all of them to be 
equal to jj,'. On the other hand we have another free parameter in the definition of negative 
powers of Pi where, for simplicity, we take j — 1/2. Since the largest negative power of Pi 
which appears in the constraint is —1/4, we will take £ — 1/2 and obtain it directly from 
Eq.(2.50), and after that we express the other negative powers by them. The eigenvalues 

for the operator are given by 

MV,Pi,P2,P3) = ^llpf^ (2.51) 

with 

h{V) = ^/V + Vc - ^/\V-V,\, and K = ^n-fXll . (2.52) 

By using these results and choosing some factor ordering, we can construct the 
total constraint operator. Note that different choices of factor ordering will yield dif- 
ferent operators, but the main results will remain almost the same. By solving the 
constraint equation Ch • ^ = 0, we can obtain the physical states and the physical 
Hilbert space "Hphys- As a final step, one would need to identify the physical observables, 
that in our case would correspond to relational observables as functions of the internal 
time (j). These steps have proven to be exceptionally difficult and have prevented from 
solving the resulting difference equations numerically, even for the simplest case of Bianchi I. 



C. Effective Equations 

When analyzing the numerical solutions of the k—0, A—0 FLRW model, the authors of 
[28] noticed that sharply peaked states followed trajectories in the (V, 0) plane that have 
a bounce, and therefore do not satisfy the classical Einstein equations. Furthermore, they 
realized that the expectation value of y|<^ does indeed follow a trajectory that satisfies (to a 
very good approximation) some equations that are now referred to as the effective equations. 
As it turns out, these effective equations can be derived from an effective Hamiltonian 
constraint Ccs- The question that arises then is how to derive, from the quantum theory 
defined by a quantum constraint C, the effective Hamiltonian. A second question pertains 
to the domain of validity of these effective equations. That is, for which states and in which 
regimes are these equations a good approximation to the exact quantum dynamics? As we 
shall see in this part, for the models that are well understood, effective equations describe 
very accurately the dynamics for appropriately defined semiclassical states. 

In the case of models for which we do not posses the full quantum dynamics, one can 
expect that the effective theory to describe very well the quantum theory for semiclassical 
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states far from the 'deep quantum regime' (where it is expected to fail). Thus, in the 
anisotropic Bianchi I, II and IX models, the effective description that we shall here consider 
provide a description in which the singularity is also replaced by a bounce. 

Let us begin by briefly describing how one obtains this effective descriptions from the 
quantum theory. The idea is to employ the geometric formulation of quantum mechanics 
[51], which provides an appropriate formalism from which one can find the effective Hamil- 
tonian constraint Ccs by computing the expectation value (C)^ of the quantum Hamiltonian 
constraint on an appropriately defined semi classical state ip. From that expression one can 
find the effective equations of motion by replacing Ceff in Hamilton's equations: q — {q, Ces} 
and p = {p,Ccs}- 

Let us now be more precise. In the geometric formulation of quantum mechanics the 
space of quantum states is seen as a symplectic space Fg, equipped with a symplectic 
structure Qq that is given by the imaginary part of the Hermitian inner product (•, •) on 
H. For each observable F one can define a function F :— (F) on normalized states. There 
is a corresponding Hamiltonian vector field for each function Xp = Q^^djiF. There is an 
interesting interplay between these vectors and the vector one would obtain by acting with 
the operator F on a state ^, 

= inXf\^ (2.53) 

Furthermore, the commutator of observables in the Hilbert space and the corresponding 
quantum Poisson bracket {F, G}q := daFdpG satisfy the relation, 

{[FA)-if^{F,G}Q (2.54) 

Thus, quantum dynamics can just be seen as ordinary Hamiltonian dynamics on the quantum 
phase space Fg, as defined by the corresponding vector field X*^. How can we relate then this 

quantum evolution with the classical evolution on the phase space F? The idea is to project 
the dynamics on Fq to F by means of appropriate coordinate functions. To be precise, let 
us assume that the classical phase space F has coordinates {q\Pi). In the Hilbert space one 
has the corresponding operators (q^^Pj). Then, one can define the projection H : Fq — > F as 
follows: H : ^ — >■ {q^,pj). One can now, given a quantum dynamical trajectory Tj on Fg, 
define the corresponding projected classical trajectory '-ft in F as: 7^ = n(T(). The question 
that arises then is whether one can find an effective Hamiltonian i^eff , defined on the classical 
phase space (and therefore being a function of {q\Pj) and possible some parameters), such 
that the trajectory 7^ = {q^,Pj) follows Hamilton's equations = {g*,ifefr} and pj = 
{pj.Hcs}- For this conditions to be satisfied, one must choose a particular 'initial state' in 
order to select a preferred trajectory T^. In practice one looks for something simpler. In 
the so called, 'embedding approach', one seeks an embedding F ^ f g C Fq of the finite 
dimensional phase space into the infinite dimensional quantum space Fq that is well suited to 
capture the quantum dynamics, in the sense that the dynamical evolution lies approximately 
within Fq. To define Fq, for any given point 7" G F, where 7° = {qi,Pi), one prescribes 
a quantum state \t'^o for all 7° G F. A first requirement is that the embedding should be 
such that q° = {^^o q^ ^/^o) and p° = {^^Pi ^7°)- second condition is dynamical and 
non-trivial: it requires that the quantum Hamiltonian vector field should be approximately 
tangent to Fq. If this is satisfied, one can project the exact quantum evolution Tt to Fq 
to obtain Tj, and from this, project down to jt — H(Tt). It is natural to regard, as a 
candidate for ifeff, the expectation value of the quantum Hamiltonian on the embedded 
submanifold: Hes{q°,Pj) :— {^joH'^^o). One should note that for the ordinary harmonic 
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oscillator, coherent states represent an exact dynamical embedding. That is, the exact 
quantum evolution lies within Tq and the effective Hamiltonian coincides with the classical 
one. There are no quantum corrections to the dynamics from these states. 

Let us now consider some important cases in homogeneous loop quantum cosmology. 



1. k=0 FLRW cosmology 

Using the geometric methods of quantum mechanics just described, one can write an effec- 
tive Hamiltonian which provides an excellent approximation to the behaviour of expectation 
values of Dirac observables in the numerical simulations [52] . The effective Hamiltonian will 
in principle also have contributions from terms depending on the properties of the state such 
as its spread. Effect of these terms turns out to be negligible as displayed from the detailed 
numerical analysis [34, 38]. Thus, the effective Hamiltonian constraint is, for A^=l, 

" 8^^ ^'"p - C^^tt , (2.55) 

which leads to modified Friedman and Raychaudhuri equations on computing the Hamilton's 
equations of motion (as we shall see below). Using (2.55) one can find that the energy density 
p = -f/^matt/^ equals 3sin^(Ab)/(87rG'7^A^). Since the latter reaches its higher possible value 
when sin^(Ab) = 1, the density has a maximum given by 

_ 3 

Thus, we see that that the maximum energy density obtained from the effective Hamiltonian 
is identical to the supremum psup for the density operator in k=0, LQC The difference is, 
of course, that in the effective dynamics every trajectory undergoes a bounce and reaches 
the maximum possible density, while in the quantum theory not every state is close to the 
critical density at the quantum bounce. 

It is easy to solve for the dynamics defined by the effective Hamiltonian. The equations of 
motion are found using the effective constraint: dtF =: F = {F,C}, with t the cosmic time. 
The only equation of motion different from the classical one (on the constraint surface) is 

V = sin (Ab) cos (Ab) , (2.57) 
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leading to the modified Friedman equation for the Hubble parameter 
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where pmax = ^ ^ is the scalar field density at the bounce. For every trajectory there are 
quantum turning points at b = where V — 0, corresponding to a bounce. Note that, at 

the bounce V\i3=^ = 2a^Kpmax > 0, so the bounce corresponds to a minimum of volume. 
Also, note that the Hubble parameter is absolutely bounded \H\ < 1/(2A7), indicating 
that the congruence of cosmological observers can never have caustics, independently of the 
matter content. 
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In the case of effective theories the proper time appears as a natural choice for an evolution 
parameter, but one can always look for internal, relational notions of time. Since, b < 
one can choose b as a relational time in the effective theories, and consider the evolution 
with respect to b. The advantage of this election is that no external time variable is needed. 
Every trajectory, that corresponds to b > 0, has a bounce at b = ^, and this value tends 
to infinity as A — 0. 

In the effective theories, we consider the interval b e [~^) should note that 

all functions and observables in Fx are periodic in b with period tt/A. It is then completely 
equivalent to regard the coordinate b as compactified on a circle. The solutions are defined 
for every t and are given by [53], 

cotAb=^, Vx{t)^^p^ V7'A2 + 9i2, (2.59) 

and 

«<)^^. + A, + i.n|y-^^. (2.60) 

a 3^0 + Vt'A^ + 9tl 

so that (f)\(to) = 00 + A<^ and the initial condition approaches the classical one (for t — to) as 
A ^ 0. Note that 0a(O) ^ ^InA as A ^ 0. Let us now consider an intrinsic description 
of the dynamics in terms of the scalar field. One can solve y as a function of 0, 

Vx{(f>) = 1/j_e"(^Snb)((/>-(;!.(to)) _|_ Y_^-a{sgab){(t>-4>{to)) ^ (2.61) 

where V+ = ^(K) + \A^F^) and V_ = ^{V+)-\ where Vo = V{(/){to)), and (3 = l^fXap^. 
Note that the effective theory recovers the quantum dynamics of {V)\^ exactly, for all states 
of the physical Hilbert space. That is, there are no further quantum corrections to the 
dynamics of Va(0). 

With this, one can see that the effective theory defines an effective homogeneous and 
isotropic spacetime metric, that takes the form, 

(ds2)eff = -dt + a^(i)efr dx^ (2.62) 

1 

with a(t)efr = (f ) ^ ^ (7^A^ + 9 1^) e . It is trivial to see that in the A — >■ limit, one recovers 
the classical spacetime metric satisfying Einstein equations. 

As we have seen, the quantum corrections captured by the effective Hamiltonian modify 
the Friedman equation in a non-trivial way, ensuring that quantum effects become important 
near the Planck scale in such a way that a repulsive force is capable of stopping the collapsing 
universe and turn it around into an expanding phase. Let us explore a little bit more how 
this quantum repulsive force can be seen. First, a modified Raychaudhuri equation can be 
written [54], 

' = -^pfl-4^V4vrGpfl-2^') . (2.63) 

It is also illustrative to write an equation for the rate of change of the Hubble parameter 
[55], 

H = -inG{p + P) ( 1 - 2 — j . (2.64) 
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These equations imply that the matter conservation equation 

p + 3H{p + P)=0, (2.65) 

has the same form as in the classical theory, even when both Friedman and Raychaudhuri 
equations suffer loop quantum corrections. From Eq. (2.64) one sees that, for matter sat- 
isfying the WEC, there is a super-inflationary phase, corresponding to if > 0, whenever 
the matter density satisfies p > Pmax/2. Note that in the A — >■ limit, we recover the 
corresponding classical equations. 

Another system of interest, for the remainder sections of this Chapter, is a scalar field 
subject to a potential V(0). Even for the simplest potential V(0) = m^</)^/2 the classical 
dynamics is drastically modified; after the big bang there is a, 'slow roll', infiationary period. 
A pressing question is how this dynamics gets modified in the effective LQC scenario. We 
know that every trajectory follows the effective Friedman equation (2.58) and has a bounce 
when p = Pmax, followed by a period of superinffation. How does that behavior affect the 
presumed inflationary period occurring at much smaller densities? First note that in that 
case, the energy density has the form: p = 0^/2 -|- m^0^/2, so there is a convenient way of 
depicting the bounce as the curve, in the (0, 0) plane, satisfying Pmax = 0^/2 + m^0^/2. The 
dynamics is therefore bounded by such ellipsoid. The equation satisfied by the scalar field 
has the same form as in the classical case: -I- SH + V^<^ = 0. One can solve these equations 
numerically [56, 57] and finds that after the superinflationary phase, the dynamics follows 
very closely the GR dynamics and exhibits an 'attractor' behavior as well. As we shall see 
in later sections, this feature of the dynamics is responsible for phenomenologically relevant 
inflation to be generic. 

Let us end this part with some comments, i) This set of effective equations have the 
property that one recovers General Relativity in the small density 'IR' limit, and that they 
are independent of the flducial V. These are non-trivial requirements that impose strong 
conditions on the particular form of the quantum constraint operator [58] . ii) Inverse volume 
effects can introduce modiflcations to the effective equations that have various consequences, 
such as loss of the universal conservation equation for matter, and extra superinffationary 
corrections. However, the physical validity of considering such inverse correction for the 
k=0 is seriously challenged, iii) It has been shown that for generic matter content, the LQC 
effective equations imply that strong singularities are generically resolved [54] . iv) A consis- 
tency check for the validity of effective equations pertains to the behavior of appropriately 
defined semiclassical states. Such states have been constructed and the predictions of the 
effective theory put to the test [37]. It was shown that both the density at the bounce and 
the minimum value of volume are very well described by the effective theory. 

2. k=l FLRW 

Let us now start with the isotropic closed FLRW model. As discussed before, there are 
two quantization available for this model. Correspondingly, the effective equations will yield 
two inequivalent theories. For the first quantization, based in the curvature as defined by 
closed holonomies, and neglecting the so called inverse triad corrections, one can arrive at 
the form of the effective Hamiltonian constraint. 



Ceff = -^-^-^V [sin\\h -D)- sin^ L> + (1 + t^)^^] + pV (2.66) 
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with D :— X-d/V^/^. We can now compute the equations of motion from the effective 
Hamiltonian as, 

V = {V^,Ceff} = {VM^ = ^^^si^l-^b - D) cos(Ab - D) . 
Prom here, we can find the expansion as, 

^ = ^ = A sin(Ab - D) cos(Ab - D) = ^ sin 2(Ab - D) . (2.67) 

^From the above equation we can see that the Hubble parameter is also absolutely bounded 
by \H\ = \0\/3 < I/2A7. We can now compute the modified, effective Friedman equation, 
by computing 

H' = ^( + sin^ D-{1 + f ) ( 1 - _ sin^ D + {1 + j')D' 
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(2.68) 

where pi = pmax[(l + 7^)-D^ — sin^ D] and pmax = 3/{87rG^'^X'^) is the critical density of the 
A; = FLRW model. 

Let us now consider the other quantization, based on defining the connection using 
holonomies along open paths. As mentioned before, this is the only available route for 
anisotropic cosmologies when there is intrinsic curvature (such as in Bianchi II and IX) . The 
effective Hamiltonian constraint one obtains from that quantum theory [40], when neglecting 
inverse scale factor effects (as was done in [38] and [60]), is 

Ceff = - 8^g^2p ^ [(^i^-^b - Df + ^'D'] + pV . (2.69) 

It is then straightforward to compute the corresponding effective equations of motion. In 
particular, by computing V — {y,Cefr}, we can find the expression for the expansion as 

^ = A cos Ab (sin Xh - D) . (2.70) 
A7 

Note that in this case, the expansion (and Hubble) is not absolutely bounded, due to the 
presence of the term linear in D. An important feature of these effective equations is that 
they describe with great accuracy the expectation value of volume during the numerical 
evolution of semiclassical quantum states [38]. It is also worth notice that for large values 
of the recollapsc volume, the effective and the classical equations coincide. In the case of 
the connection based quantization [40], there are two different bounces, that approach the 
unique bounce of the curvature based equations when the universe grows to be large [40]. 
Let us now consider the effective equations for anisotropic models. 
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3. Anisotropic Models: Bianchi I, II and IX 



Considering the effective description of anisotropic models is interesting in view of tlie 
BKL conjecture [61, 62], tliat states tfiat locally, generic spacetimes approaching the clas- 
sical singularity behave as a combination of Bianchi cosmological models. The effective 
Hamiltonian constraint for Bianchi I and II can be written in a single expression [47, 49, 63] , 
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where the parameter e allows us to distinguish between Bianchi I {£= 0) and Bianchi II (e— 
1). This Hamiltonian together with the Poisson Brackets {d,Pj} — SnG^Sj and {(f), p^p} — 1 
gives the effective equations of motion. In these previous effective Hamiltonians, one chooses 
the lapse N = V. 

In Bianchi IX, we choose A^=l to include more inverse triad corrections, then the effective 
Hamiltonian is given by [63] 
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Let us discuss the issue of singularity resolution when these equations are studied numeri- 
cally, i) All solutions have a bounce. In other words, singularities are resolved. In the closed 
FRW and the Bianchi IX model, there are infinite number of bounces and recollapses due to 
the compactness of the spatial manifold, ii) One can have a different kind of bounce dom- 
inated by shear cr, but only in Bianchi II and IX. In Bianchi I, the dynamical contribution 
from matter is always bigger than the one from the shear, even in the solution which reaches 
the maximal shear at the bounce [66]. iii) In the flat isotropic model all the solutions to 
the effective equations have a maximal density equal to the critical density, and a maximal 
expansion (6'^a,x = STrGpmax = 3/(27A)) when p = Pcrit/2. For FRW k = I model, every 
solution has its maximum density but in general the density is not absolutely bounded. In 
the effective theory which comes from connection based quantization, expansion can tend 
to infinity. For the other case, expansion has the same bound as the fiat FRW model. How- 
ever, by adding some more corrections from inverse triad term, one can show that actually 
in both effective theories the density and the expansion have finite values, iv) For Bianchi I, 
in all the solutions p and 9 are upperly bounded by its values in the isotropic case and a is 
bounded by cr^ax = 10.125/(37^A^) [64, 65]. For Bianchi II, 9, a and p are also bounded, but 
for larger values than the ones in Bianchi I, i.e., there arc solutions where the matter density 
is larger than the critical density. With point-like and cigar- like classical singularities [66], 
the density can achieve the maximal value {p ~ 0.54ppi) as a consequence of the shear being 
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zero at the bounce and curvature different from zero, v) For Bianchi IX tlie befiaviour is 
tlie same as in closed FRW, if tlie inverse triad corrections are not used, then the geometric 
scalars are not absolutely bounded. But if the inverse triad corrections are used then, on 
each solution, the geometric scalars are bounded but there is not an absolute bound for all 
the solutions [63, 64]. 

III. INHOMOGENEOUS PERTURBATIONS IN LQC 

The theory of quantized fields in curved space-times has become an essential tool in mod- 
ern early-universe cosmology. In that framework one studies the behavior of quantum fields 
propagating in space-times with generic Lorentzian geometries, as in General Relativity. 
One expects this theory to describe accurately physical processes in situations where we are 
confident about the validity of its building blocks: a description of matter fields in terms of 
quantum field theories, and a space-time geometry given by a smooth, classical space-time 
metric. These assumptions are reasonable, for instance, during the inflationary era in which 
the energy density and curvature are believed to be more than ten orders of magnitude 
below the Planck scale. However, earlier in the history of the universe, closer to the Planck 
era, quantum gravity effects become important and the description of space-time geometry 
in terms of a smooth metric is expected to fail. To include physics in the Planck regime 
QFT in curved backgrounds needs to be generahzed to a QFT in quantum space-times. The 
singularity-free quantum geometry provided by LQC, summarized in the previous section, 
provides a suitable arena to formulate such a theory, and the quantization of scalar fields 
on those quantum cosmologies was introduced by Ashtekar, Kaminski and Lewandowski in 
[3], and further developed in [4, 67-69]. Having in mind the most interesting application of 
this framework, we summarize here the construction of the QFT of scalar and tensor metric 
perturbations propagating in a quantum FLRW universe, i.e. the quantum gravity theory of 
cosmological perturbations. For more detailed information, see [3, 4]. 

As mentioned in the introduction of this chapter, the construction will follow the guiding 
principle that has been useful in the quantization of the background: first carry out a trunca- 
tion of the classical theory to select the sector of General Relativity of interest, and then move 
to the quantum theory by using LQG techniques. Starting from General Relativity with a 
scalar field as matter source, we will truncate the phase space to the sector containing cos- 
mological backgrounds plus inhomogeneous, gauge invariant, first order perturbations, and 
then write down the dynamical equations on that classical, reduced phase space. The main 
approximation behind this truncation, and underlaying the subsequent quantization, is that 
the back-reaction of inhomogeneous perturbations on the homogeneous degrees of freedom 
is neglected. The second step is to move to the quantum theory. Physical states will depend 
on background homogeneous degrees of freedom as well as on inhomogeneous ones. Our 
basic approximation, however, enables us to write these quantum states as tensor product 
of the homogeneous part, which will evolve independently of perturbations, and first order 
inhomogeneities thereon. The homogeneous part will therefore be the same as the quantum 
geometries obtained in the previous section, in which the big bang singularity is replaced by 
a bounce. The surprising result appears in the evolution of perturbations. Without further 
approximation, the evolution of inhomogeneities on those quantum geometries turns out to 
be mathematically equivalent to the quantum theory of those fields propagating on a smooth 
background characterized by a metric tensor. The components of that smooth metric, how- 
ever, do not satisfy the classical Einstein equation. They are obtained from expectation 
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values of certain combinations of background operators, and incorporate all the informa- 
tion of the underlying quantum geometry that is 'seen' by perturbations. The message is 
that the propagation of inhomogeneous perturbations is not sensitive to all the details of 
the quantum space-time, but only to certain aspects, which appear precisely in a way that 
allows to encode them in a smooth background metric. This is an unforeseen simplification 
that facilitates enormously the treatment of field theoretical issues. 

The last step is to develop the necessary tools to check the self-consistency of this 
construction. It is necessary to show that, in the physical situations under consideration, 
the Hilbert space of physical interest contains a large enough subspace in which the 
back-reaction of perturbations on the background is indeed negligible, in such a way that 
our initial truncation is justified. That should be done by comparing the expectation value 
of the Hamiltonian and stress-energy tensor for perturbations with that of background 
fields. Those computation will require of techniques of regularization and renormalization. 



A. The classical framework 

The goal of this subsection is to summarize the construction of the truncated theory 
of classical FLRW space-times coupled to a scalar field, plus gauge invariant, linear per- 
turbations on it, and write down the equations describing their dynamics. The reader is 
referred to the extensive literature for more details (see, for instance, [70]). Wc adopt here 
the Hamiltonian framework which, as shown in [71], is particularly transparent on the task 
of finding gauge invariant variables. It will also provide the appropriate arena to pass to the 
quantum theory in the next section. For simplicity and for physical interest, we work here 
with a spatially flat FLRW universe. 

The procedure can be divided in three steps: 1) Starting from the full phase space, expand 
the configuration variables and their conjugate momenta in perturbations, and truncate the 
expansion at first order. Expand also the constraints of the theory (the scalar and vector 
constraints) and keep only terms containing zero and first order perturbations. 2) Use 
the constraints linear in first order perturbations to find gauge invariant variables. Those 
variables coordinatize the so-called truncated, reduced phase space. 3) Use the part of the 
constraints quadratic in zero and first order perturbations to write down the dynamics. See 
[4] for further details and subtle points of this construction. 

1. The truncated phase space 

Let us consider General Relativity coupled to a scalar field on a space-time manifold 
M = S X R, with S = R"^. Due to the infinite volume in S, spatial integrals of homogeneous 
quantities will introduce infrared divergences. To be able to write meaningful mathematical 
expression, it is convenient to introduce a fiducial cell V and restrict all integrals to it. V 
can be chosen to be arbitrarily large, or at least larger than the observable universe. At 
the quantum level this will be equivalent to restrict to V the support of test functions in 
operator valued distributions. 

We will work with ADM variables for the gravitational sector, where the canonical con- 
jugated pairs consist in a positive definite 3-metric on E, Qab, and its conjugate momen- 
tum p"'' (the same analysis can be done in connexion variables, by including the corre- 
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spending Gauss constraint; see [72-75], [4]). The full phase space F consists of quadruples 
{qab{x),p"'^{x),^{x),Il{x)} G r, where 11(5') is the conjugate momentum of the scalar field 
Because we are interested in expanding around Fhom C F, the (FLRW) isotropic and 
homogenous sector of F, it is convenient to introduce a fiducial flat metric Qab, and use it to 
raise and lower indices. We will denote x — {xi, X2, x^) the Cartesian coordinates defined by 
Qab onV,V the volume of V with respect to Qab, which we take equal to one to simplify the 
notation, and q = 1 the determinant of Qab- 

Consider now curves 7[e] in F, which pass through Fhom at e = 0. Expanding the phase 
space variables around e = 0, we have: 

qab[e]{x) = a^Qab + e ^(labi.^) + ■■■ + —, Sq^^\x) + ... 
D a 

$[e](x) = (p + € ip^^\x) + . . . , 

U[e]{x) = p(<^) + e7rW(f) + ... (3.1) 

It is convenient to consider the first order perturbations dq''^fj{x), 6p"'^^^\x),Lp^^\x),7i^^\x) 
as purely inhomogeneous functions of x, in the sense that the integral of any of them on V is 
zero. By truncating the above expansions at first order we obtain the truncated phase space, 
made of four pairs of conjugate variables: FTmn = {(a, Pa, 4'iP{4>)i ^Qab, 5p°'''^^\ (/?^^\ tt^-*-^)} = 
rhom X Ti, where the only non-zero Poisson brackets between the basic variables are: 

{a, Pa} = 1, {Sq^^ix,), 5p'^^'\x2)} = 51 5t) 5(f i, x^), 

{<P, = 1, W^'\x,), 7r«(f2)} = 5{x^,X2), (3.2) 

where 5{xi,X2) — {5{xi,X2) — 1) is the Dirac delta distribution on the space of purely 
inhomogeneous fields. /,From now on we will work only with first order perturbations, so we 
will omit the super-index (1) to simplify the notation. 

Because of the homogeneity of the background it is convenient to Fourier transform the 
perturbation fields and carry out the standard scalar-vector-tensor decomposition, in which 
the 6 degrees of freedom of Sqab are decompose into two scalar, two vector, and two tensor 
modes (see e.g. [71], [4] for details). Because perturbations are inhomogeneous, the restriction 
to the fiducial cell V is not strictly necessary, and one can avoid the artificial quantization 
of k that it introduces. However, from the physical point of view one can absorb modes 
with wavelength larger than the observable universe in the background. Therefore, we will 
consider that the Fourier integrals incorporate an infrared cut-off ko provided by the size of 
the observable universe. 



2. Constraints and reduced phase space 

A similar expansion to (3.1) can be carried out for the constraints. In General Relativity 
the Hamiltonian is a sum of constrains, the famihar scalar S[A'"], and vector Y[N] constraints. 
If 7[e] is now a curve that lies in the constraint hyersurface of F, and intersects Fhom at e = 0, 
by referring to the constraints collectively as C{q"''',Pab, II) (suppressing the smearing fields 
for simplicity), we expand around e = to obtain a hierarchy of equations: 

CW:=CUo = 0, CW:=^U = 0, ... C^") := = 0, ... (3.3) 
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• The zeroth-order constraint, = 0, is just the restriction of the full constraint to the 
homogeneous subspace Fhom- The zeroth-order vector constraint is trivially satisfied 
because of the gauge fixing on the zeroth-order variables, introduced by the use of 
the fiducial metric qab in (3.1). The zeroth-order scalar constraint. So, is quadratic in 
zeroth-order variables and can be interpreted as the generator of background dynamics. 
This dynamics is exactly the same as that of the unperturbed theory. 

• First order constraints are linear in first order variables. They generate gauge trans- 
formations in FTrun and, as usual, tell us that some of our degrees of freedom are not 
physical. Initially we have 6 (xoo) degrees of freedoms in 6qabix), plus 1 degree of free- 
dom in the scalar field (p{x), a total of 7. As mentioned above, Sqab{x) is conveniently 
decomposed in Fourier space into two scalars, two vector, and two tensor modes. We 
have the scalar and three vector constraints, a total of 4. Therefore, the number of 
physical degrees of freedom is 7 — 4 = 3. There is an elegant systematic procedure 
to construct gauge invariant variables out of those 3 degrees of freedom, and we refer 
the reader to [71] for details. It can be summarize as follows. In FLRW backgrounds, 
scalar perturbations are affected by the scalar constraint and only one of the vector 
constraints; they reduce the three scalar degrees of freedom that we have initially, 
two from gravity and one from the matter sector, to only one physical scalar mode. 
Vector perturbations are affected by two of the vector constraints that kill completely 
the vector modes. In other words, in the absence of matter with vector degrees of 
freedom, as in the case we arc studying, there arc no physical vector perturbations. 
Tensor modes are not affected by any of the constraints and therefore the two original 
tensor modes are the physical ones, i.e. they are gauge invariant. In summary, after 
imposing the constraints we are left with one scalar degree of freedom, which we choose 
to be the familiar Mukhanov variable Q, and two tensor modes T*^^^ and T^^^. They 
are gauge invariant variables, and together with their conjugate momenta form the 
reduced, truncated phase space of first-order perturbations, FTrun- Equations C*-"^ — 
with n > 1 do not add further constraints on first oder perturbations. 

• The second-order constraints in the full phase space F involve terms quadratic in first- 
order perturbations as well as linear terms in second-order perturbations. When a 
second order constraint C^^^ is restricted to the truncated phase space FTrun, terms 
containing second order perturbations are disregarded, and the resulting combination 

of quadratic terms in first-order perturbation with coefficients containing background 
quantities, C-^-*, is no longer a constraint. The truncated second-order scalar constraint 
§2 is interpreted as the Hamiltonian that generates the dynamics of gauge invariant 

first-order perturbations. It has the form S2 = ^2 -l- §2 + ^2 , which indicates 
that scalar and tensor modes evolve independently of each other, where 

with K, = SttG. The two tensor modes behave identically, and we have denoted them 
collectively by T. For pedagogical reasons we only write down the expressions for 
tensor perturbations. See [4, 76] for explicit expressions for scalar modes. In the above 
equations are the conjugate momenta of 7g, with Poisson brackets {7^,P^,} = 
(27r)^5(^ — k'). Tensor perturbations, except for the constant factor l/(2y^) that 
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provides the appropriate dimensions, behave exactly as masslcss, free scalar fields 
(scalar perturbations behave as a scalar field subject to a time dependent 'emergent' 
potential). The (homogeneous) lapse function indicates the time coordinate one is 
using. For instance, N — 1 corresponds to standard cosmic time t, N — ato conformal 
time T), and N = a? to choosing the scalar field as a time variable, which turns 
out to be the natural choice in the quantum theory. 

To summarize, the phase space of physical interest is the reduced, truncated phase space 

f made of elements {(a, P„, 0,p(^)); (Q,^, pf\ rS'\pf''\ ri'\ pf '''))} G f ^.n- The ho- 
mogenous degrees of freedom evolve with the zeroth-order Hamntonian. This evolution takes 
place entirely in Fhom, and is independent of perturbations, refiecting the main approxima- 
tion of the truncated theory. The homogenous dynamical trajectory can then be 'lifted' 
to FTrun, providiug a well-defined evolution of first-order perturbations on the homogenous 
background. This evolution is specified by the Hamiltonian §2- 

B. Quantum theory of cosmological perturbations on a quantum FLRW 

1. Quantization o/FTrunc 

In this section we pass to the quantum theory starting from the reduced, truncated phase 
space Fttuii- The structure of the classical phase space Fttuii = ^hom x suggests that in 
the quantum theory the total wave function ^ has the form 

^(a, Tg, 0) = ^o{a, 0) ® ^(Tg, 0) . (3.5) 

This product structure is maintained as long as the test field approximation holds. Because 
back-reaction is neglected, the background part evolves independently of perturbations, 
and the solutions for are the ones obtained in section II. When written in terms of 
the relational time 0, they satisfy the equation P(^)^o = —ifid^'^Q = Ho'^o, where the 
operator = a/© is obtained from expressions (2.19) and (2.21). The remaining task 
is to 'lift' this trajectory to the full Hilbert space, by writing down the quantum theory 
for propagating on the quantum geometry specified by "^q. The evolution of ip will be 
specified by the operator analogue of §2 ■> which generates the dynamics in the classical 
phase space. In the classical theory depends on inhomogeneous degrees of freedoms, 
but also on the homogeneous ones via the scale factor a. Therefore, in the quantum theory 
the corresponding operator will act on perturbations ip as well as on ^o- 

Our goal is to generalize the theory of QFT in curved space-times in which, on the one 
hand, quantum fields propagate in an evolving classical FLRW specified by ac\{r)) and, on the 
other hand, perturbations are commonly quantized using the Heisenberg picture. Therefore, 
to facilitate the comparison, we pass in this section to the Heisenberg picture. In obtaining 
the evolution equations for the operator 7^ and its conjugated momentum we will use as 
internal time, because it is the evolution variable that appears naturally in the quantum 
theory, while standard cosmic or conformal time are represented by operators. Internal time 
corresponds to use the lapse function = a'^/p(^) in the expression (3.4). By choosing 
an appropriate factor ordering to convert it to an operator, we have (as it is common in 
quantum theory, we are not free of factor ordering ambiguities) 
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d,pri<P) = ^ipPM^] = -^(P(,f a^(0)p-,f 0r^). (3.6) 

These equations involve background operators as well as perturbations. However, the test 
field approximation allows us to 'trace over' the background degrees of freedom. This can 
be done by taking expectation value with respect to the background wave function (in 
the Heisenberg picture) obtained in the previous section 



k ' 

where background operators have been replaced by expectation values and, additionally, 
we have used the evolution equation P((/,)\E'o = -f^o^o- The test field approximation ensures 
that we are not losing any information when passing from (3.6) to (3.7). These are the 
Heisenberg equations for perturbations, in which the coefficients are given by expectation 
values of background operators in the quantum geometry specified by \I'o. This is a quantum 
field theory of cosmological perturbation on a quantum FLRW universe. Note that the 
above equation is exact, and not further approximation has been made beyond the test field 
approximation. 

In this theory, space-time geometry is no described by a unique classical metric, it is 
rather characterized by a probability distribution ^'o that contains the unavoidable quantum 
fiuctuations. The propagation of perturbations is sensitive to those fiuctuations, and not only 

the mean effective trajectory (a). However, it is remarkable that those effects can be encoded 

^ _i ^ _i 

in a couple of expectation values of background operators: {H~^) and {Ho ^ ci^{(t>) Ho ^) [3, 4]. 
Borrowing the analogy from [4], this is similar to what happens in the propagation of light 
in a medium: the electromagnetic waves interact in a complex way with the atoms in the 
medium, but the net effect of those interactions can be codified in a few parameters, such 
as the refractive index. Similarly, although the final equations (3.7) depend in a simple way 
on the quantum geometry, it had be very difficult to guess the precise 'moments' of the 
quantum geometry that are involved in the evolution of perturbations. 

We can now compare the above evolution equations with the familiar quantum field 
theory of cosmological perturbations on classical FLRW geometries, in which the Heisenberg 
equations, when is used as time, are written in terms of the classical background quantities 
a{(f)) and p(^) as 

Comparing with (3.7) we see that the QFT in a quantum background is indistinguishable 
from a QFT on a smooth FLRW metric 

gabdx'^dx^ = ds^ = -(p(^)) -2 5*^(0) d(j)^ + a{(j)fdx^ (3.9) 
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where 

(m) = (^o ) and a = . (3.10) 

K^o I 

In terms of the more famihar conformal time used in cosmology, we have ds^ = a^{fi) (— d?7^ + 
da^^), with dfj = [a^(0)] p^^^ dep. This smooth metric captures all the information of quantum 
geometry that is 'seen' by perturbations. Note that its components contain h and it does 
not satisfy the Einstein equation, not even the LQC effective equations. 

In terms of this smooth metric, we can write the Heisenberg equations (3.7) as a second 
order differential equation 

'H' + ^~^rS + k^%^0, (3.11) 

where the prime now denotes derivative with respect to fj. This equation is mathematically 

equivalent to the familiar formulation of QFT in classical FLRW space-time, where all the 
effects of the quantum background geometry have been encoded in a dressed, smooth met- 
ric tensor Qab- This unexpected mathematical analogy highly simplifies the analysis, not 
only conceptually, but also at the technical level. It allows to extend well-stablished tech- 
niques from classical space-times to define the physical Hilbert space and the appropriate 
regularization and renormalization of composite operators on it (see [4, 76] for details of 
that construction). These are the necessary tools to make sense of the momentum inte- 
grals appearing in, e.g. the Hamiltonian S2, that so far were formal, and to regularize the 
expectation value of the energy-momentum tensor in the physical Hilbert space. 



2. The physical Hilbert space 

In this subsection we briefly summarize how techniques of regularization and renormaliza- 
tion from linear QFT on classical space-times, can be extended to characterize the physical 
Hilbert space of cosmological perturbations on quantum backgrounds, and to regularize 
composite operators on it. Among the existing methods of regularization we will work in 
the adiabatic approach [85, 87], which is particularly convenient to perform explicit compu- 
tations, including the numerical implementation required in the next section. 

The spatial homogeneity and isotropy of our FLRW background allows us to expand the 
field operator T(a;, fj) in Fourier modes (a similar construction holds for scalar perturbations) 

f{x,fj) -j^.Jd'k (4-e.(77) +4e^(77)) e^'^ . (3.12) 

The field operator ^{x^f]) satisfies the equation of motion (3.11) as long as the mode func- 
tions ek{fj) are solution of the wave equation 

4(^) + 2?^4(^) + ^'e,(77)=0, (3.13) 

(X 

were prime indicates derivative with respect to fj. The solutions ek{fi) can be understood as 
'generalized positive frequency modes', because they play the role of standard positive fre- 
quency solutions e~*^*/ \/2k in Minkowski space-time. The canonical commutation relations 
for the field operator T(x, fj) and its conjugate momentum imply 

[4,4,] = ih{27:f5{k-k') {ek(fj),e,,(fj))-' ; [4,4] = 0, (3.14) 
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where 



{ek{fi),ey{fi)) := —{ek{f,)e'*,{f,) - e'M^'ifl)) ■ (3.15) 



Therefore, if we impose the normahzation condition {ek{fi), ek{f})) = i, and At will satisfy 
the familiar commutation relations of creation and annihilation operators. Note that the 
scalar product {ek{v),ek'{fj)) is constant in time if ek{fj) and 6^/(77) are solutions of (3.13). 
The Hilbert space is then constructed as follows. The vacuum state |0) (associate with the 
set of generalized positive frequency modes e^) is defined as the state annihilated by all 
Aj^. The associated Fock space "Hi arises by the repeated action of creation operators At 
on the vacuum. It is important to notice that the vacuum state constructed in this way is 
translational and rotational invariant, as can be checked, e.g. by explicit construction of the 
two point function. 

It is clear from the construction that a different choice for the generalized positive fre- 
quency bases in (3.12), provides different and At operators, and therefore a different 
definition of vacuum state. None of those vacua is preferred as compare to the others. Even 
more, different vacua may not even belong to the same Hilbert space, and the quantum 
theories constructed from each of them arc in that case unitarily incquivalcnt. 

The existence of unitarily inequivalent quantization is common in QFT in curved space- 
times (see e.g. [84]). In cosmological backgrounds, however, it is possible to add appropriate 
regularity conditions to the mode functions to select a preferred Hilbert space. The 
adiabatic condition [85-87] in FLRW backgrounds imposes that, in the asymptotic limit in 
which the physical momentum k/a is much larger than the energy scale provided by the 
space-time curvature Er, must approach the Minkowski space-time positive frequency 
modes, e~^^^ / ^/Tk: . at an appropriate rate (for a brief summary see, e.g. [4], and references 
cited there). The modes satisfying this conditions are called modes of Nth adiabatic 
order, and the associated vacuum an adiabatic vacuum of the same order, where the order 
is specified by the exact rate of approach to the Minkoswkian solutions. 

Notice that the adiabatic condition does not single out a preferred vacuum, because there 
are many different families satisfying it to a given order (it imposes only an asymptotic 
restriction for large {k/a}/ Er). However, it is possible to show that if we restrict to adiabatic 
order N > 2, all different vacua belong to the same Hilbert space "Hi. (This is strictly true 
if we restrict our QFT to the compact fiducial cell V. In the non-compact case one needs 
to be more precise in the sense in which the Hilbert space is unique [84], because infra-red 
divergences appear). Additionally, if > 4 there is a well defined procedure to extract the 
physical, finite information from the formal expression of operators of interest for us, the 
Hamiltonian and the stress-energy tensor, by subtracting ultra-violet divergences in a local 
and state independent way, while respecting the covariance of the theory. 

The Hamiltonian operator generating time evolution (in conformal time) , and the energy 
density p, are related by 



2,formal (27r)^ 



If |0) is a 4th-ordcr adiabatic vacuum associated with a family of solutions 6^(77), the renor- 
malized expectation value of the energy density is given by 



8^527 (27r) 



(3.17) 
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with 

= _ + _ , (3.18) 

where C{k,ri) are the subtraction terms provided by adiabatic regularization [87]. The 
renormahzed expression for the expectation value of the Hamiltonian operator is obtained 
from the previous three equations. The above subtractions make the expectation value of 
the hamiltonian and energy density finite for any state in the Hilbert space of 4th-order 
adiabatic states, "Hi. Additionally, the procedure has the properties that any method of 
regularization/renormalization are expected to satisfy, enunciated in the Wald's axioms [84]. 
Although strictly speaking the above expressions provide only quadratic forms in the Hilbert 
space, recent results indicate that they are expectation values of operator value distributions 
pC^) and S^P in Tii. 

In summary, our QFT in quantum FLRW admits a straightforward extension of the 
adiabatic approach of hnear QFT in classical backgrounds. The physical Hilbert space Hi 
is then singled out by restricting to 4th order adiabatic states. In addition, the adiabatic 
condition provides the necessary control on ultra-violet divergences that allows a systematic 
procedure to regularize the Hamiltonian and the stress-energy tensor on "Hi. This completes 
the formulation of the theory. 

3. Criterion for self- consistency 

The last step in the construction is to check whether the underlaying approximation in 
our truncated theory, the test field approximation, is satisfied throughout the evolution. In 
our QFT in quantum space-times this question translates to check whether the expectation 
value of the stress-energy tensor can be neglected when compared to the background one. 
However, in an homogeneous and isotropic background a sufficient condition for this to be 
satisfied is that energy density on scalar and tensor perturbations, (p(^)), be much smaller 
than the background energy density (po) at any time during dynamical phase of interest [4]. 
It is evident that one can always find states for perturbation for which that requirement is 
not satisfied. Therefore, the relevant question is: is there a sufficiently large subspace of the 
physical Hilbert space for which the previous condition on the energy density is satisfied? 
If the answer is in the affirmative then one has a self-consistent approach in which test- 
field approximation holds. This is a key question to ensure self-consistency, and has to be 
answered when this framework is applied to a concrete physical problem, as we do in the 
next section. 

C. Comments 

The previous framework is suitable to face interesting conceptual questions arising in 
quantum gravity. For instance, when does standard QFT in curved space-times become 
a good approximation? Is it safe to use standard QFT during inflation? This question 
can be answered straightforwardly because both theories have been written in the same 
form. From equation (3.11) it is clear that the standard QFT is recovered in the regime in 
which the quantum aspects of the geometry can be neglected, and Section (II C) provided 
the conditions under which this happens. When the background energy density (po) is 
below one thousandth oi pp£, quantum corrections become negligible and General Relativity 
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becomes an excellent approximation. This is the regime in which standard QFT arises from 
the more fundamental framework presented in this section. Therefore, in the inflationary 
era where (po) < lO^^^pp^, we expect the familiar QFT to be an excellent approximation. 

By construction, this framework encompasses the Planck regime and is suitable to discuss 
trans-Planckian issues and distinguish real problems from apparent ones. In LQG there is 
a priori no impediment for trans-Planckian modes to exist. It may seem at first that the 
existence of a minimum area may preclude their existence, but quantum geometry is subtle 
and, for instance, there is no minimum value for volume or length. In addition, if we pay 
attention to the construction of the background quantum theory, trans-Planckian quantities 
appear there without causing problems: the value of the momentum p(^) of the background 
scalar field (j) is generally large in Planck units. However, the background energy density 
is hounded above by a fraction of the Planck energy density. Something similar happens in 
our quantum field theory. There trans-Planckian modes are admitted as long as the total 
energy density in perturbations remains small as compared to the background. That is the 
real trans-Planckian problem, which becomes a non-trivial issue in the deep Planck regime 
where the volume of the universe acquires its minimum value. 

IV. LQC EXTENSION OF THE INFLATIONARY SCENARIO 

The previous sections have summarized the physical ideas and mathematical tools nec- 
essary to undertake the quantization of the sector of General Relativity containing the 
symmetries of cosmological space-times and the study of cosmic perturbation thereon. The 
goal of this section is to apply those techniques to extend the current picture of the evolution 
of our universe to include the Planck regime. 

The cosmological ACDM model with an early phase of infiation contains conceptual lim- 
itations that are dictated by the domain of applicability of the physical theories in which it 
is based: General Relativity and Quantum Field Theory. One needs a theory of quantum 
gravity to extend the model to include physics at the Planck era. Subsection IV A summa- 
rizes how, by introducing a scalar field with suitable potential, LQC provides a space-time 
in which the big-bang singularity is resolved by the quantum effects of gravity, and in which 
an inflationary phase arises almost unavoidably at later times. In subsection IV A it is 
shown how the evolution of cosmological perturbation can be extended to include the pre- 
infiationary space-times provided by LQC. In this sense the current scenario for the evolution 
of our universe and the genesis of cosmic inhomogeneities is extended all the way to the big 
bounce [78]. This extension goes beyond the conceptual level, as it appears a narrow window 
in which the effects of Planck scale physics could be imprinted in the CMB and galaxy dis- 
tributions, and concrete ideas connecting those effects with forthcoming observations have 
been proposed. 

A. Inflation in LQC 

As we have mentioned in previous sections, after the bounce there is a period of superin- 
fiation where H > until the density reaches half its value at the bounce, after which one 
has H < 0. It was first hoped that this period of superinfiation would be enough to account 
for the necessary number of e-foldings compatible with observations, but this period turns 
out to be too short when there is no potential for the scalar field. Thus, it is clear that 
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one needs such a potential to compare the LQC predictions with the inflationary paradigm. 
The simplest case one can consider is quadratic potential V(0) = (l/2)m^0^, that has been 
extensively studied in the literature and is compatible with the 7-years WMAP observations 
[77] . The existence of the bounce solves one of the conceptual challenges that the standard 
scenario, based on the GR dynamics poses. That is, in the GR dynamics, there is always 
a past singularity, even in the presence of eternal inflation [5]. The standard formalism is 
therefore, conceptually incomplete. 

The question that we shall pose in this part is the following: Can we estimate how 
probable it is to have enough inflation for the cosmological background? Let us be more 
precise with the question. We know that every effective trajectory undergoes a bounce, and 
some of them will experience enough e-foldings and will be of phenomenological relevance. 
Rather amazingly, WMAP has provided us with a small observational window for the scalar 
field at the onset of inflation [77, 79], written in terms of a reference time t^^ for which a 
reference mode used by WMAP exited the Hubble radius in the early universe. With an 
4.5% accuracy, the data is, in Planck units [77, 79]: 

= ±3.15 , 0(tfej = Tl-98 X 10-^ , H{tk,) = 7.83 x 10"^ . 

We can now pose the question more precisely. From all the solutions S to the effective 
equations in LQG, how many of them pass through the allowed interval? This poses yet 
another question. How are we going to 'count' trajectories? Is there a canonical way of 
measuring them? A proposal to answer this question was put forward long ago [80, 81] 
based on the idea of using the Liouville measure on phase space S, that is invariant under 
time evolution. The idea then is to compute the volume of §wmap, those solutions that pass 
through the WMAP window, relative to the total volume of S: 

P.„b=M^. (4.1) 

In order to compute this probability, one has to be careful with the way one measures all 
possible trajectories (for a discussion sec [82]). 

Let us begin with the kinematical phase space (V, b; 0,p(<^)). The constrained surface f 
(as defined by the constraint C) is three dimensional and can be given coordinates {V, b, 0). 
But in that surface, the symplectic structure is degenerate and does not define a volume 
form. For that one has to go to the space of physical states, or reduced phase space P, 
formed by the gauge orbits on the constrained surface. An alternative is to perform a gauge 
fixing to select a two dimensional surface which is transversed only once by each gauge 
orbit. As we have seen in the previous section the evolution of coordinate b is monotonous, 
so one can fix the gauge by selecting b = bo- With this choice, P has coordinates (V,0). 
Now, the pullback ^2 of ^2 to P defines the Liouville measure there. The problem is that, 
with respect to this measure, the volume of P is infinite! One has to define a procedure to 
'regularize' the integral to have finite results. The key observation is that, in the A;=0 case 
we are considering, there is an extra gauge freedom that arises from the fact that the size 
of the fiducial cell V one starts with is arbitrary. This means that a rescaling of the cell 
V — )■ i^V should leave the physics invariant. This rescaling translates into a rescaling of the 
canonical variables as V ^ i^V, (?P{<t>)-i while b and 4> remain invariant. However, 

this transformation on phase space does not leave the symplectic invariant, so it can not be 
regarded as a canonical transformation. Still, one has to gauge out this symmetry in order 



38 



to obtain truly physical quantities. The problem is that Q does not project down to the 
quotient. One possibility is to preform a further 'gauge fixing' by selecting a cross section § 
of S. For instance, one could choose a given value of volume V = Vo (for our previous choice 
b = bo). One can then restrict the measure to the cross section S to obtain the measure d/i, 
that now depends only on 0. In effective LQC, one has 



dp, 



37r 

^ sin2(Abo) - 87r^2V(0) 



1/2 

d0, (4.2) 



As expected, the GR measure is obtained by taking A — ?■ [57]. Even when the construc- 
tion involved the Liouville measure that is invariant under Hamiltonian time evolution, the 
resulting measure djl on the space of physically distinguishable configurations S depends on 
the choice of gauge fixing parameter bo, in a non-trivial way [57, 79]. A choice of b, in turn, 
fixes a value of the energy density p = po, which implies that the probability will depend on 
the energy density at which it is computed. In General Relativity there is no natural value 
of density for computing the probability, other than the big bang itself. The problem is that 
the density is infinite there and the range of is unbounded, so the volume is also infinite. 
Another possibility would be to introduce a cut-off at, say, the Planck density [83], but there 
is no reason to believe that GR is valid at that scale. In fact, one of the main lessons of loop 
quantum cosmology is that GR is not valid near the Planck scale (in energy density) but 
the isotropic degrees of freedom are rather described by the effective LQC theory. In this 
description, there is a natural preferred density which is precisely the density at the bounce 
Pmax- Thus, in what follows we shall take the bounce as the natural point where to compute 
probability. The corresponding 'gauge fixing' implementing this choice is then bo = tt/2X. 

Let us now rephrase the question that we initially posed at the beginning of this part: 
What is the relative number of solutions Swmap that pass through the observational WMAP 
window, from the total number of solutions § at the bounce? As explained before, the 
probability is computed using formula (4.1), where the volume is now obtained by integrating 
a uniform distribution (as a function of 0). The key to computing the probability is then 
a detailed knowledge of the global dynamics, for all possible values (pB of the scalar field 
at the bounce. Extensive numerical evolutions have shown that almost all trajectories fall 
within the observational window. It is only for the small window —5.46 < (pB < 0.934 from 
the total range of 0^ G [—7.44 x 10^, 7.44 x 10^] that the future dynamics lies outside the 
WMAP window [79]. For this interval, the probabihty that the dynamics falls outside of 
the observational window is less than 3 x 10~^. To understand this, one can see the LQC 
dynamics as shown in Fig. 1, where one considers a uniform distribution at the bounce and 
follows the dynamics. As can be easily seen, most trajectories funnel into a very small region 
that is precisely where the WMAP window is. Just before the onset of inflation the density 
is approximately 10"^^ smaller than the density at the bounce. At that density the allowed 
WMAP region is only 4% of the total allowed range in [79]. Thus, as seen in the Figure, 
almost all of the trajectories starting at the bounce scale fall into a very small region at the 
onset of inflation [57] . 

One should also note that this attr actor feature of the global dynamics, together with 
the non-invariance of the measure djl, explains why the probability is much smaller when 
computed in General Relativity at the onset of inflation [57, 59]. 

Let us summarize. In LQC it is natural to consider the bounce as the point where to 
compute probability of inflation. The global dynamics is such that most of the trajectories 
get funnelled into the small WMAP window at the onset of inflation where the density is 
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FIG. 1: In this figure we plot the exterior, maximal density surface pmax and a surface of constant 
density Ponset ^ Pmax (not drawn to scale, of course) on the ((/>, (p) plane. Trajectories with a 
uniform distribution at the LQC bounce ellipsoid are plotted. Note that trajectories for which 
there is enough inflation get funnelled into a small region in the smaller ponset ellipse. Near this 
surface, the GR and LQC dynamics almost coincide 

11 orders of magnitude smaller than the density at the bounce. Thus, one can conclude 
that having enough inflation is generic in loop quantum cosmology for the homogeneous and 
isotropic background, when semiclassical states are considered. 

B. Pre-inflationary evolution of cosmic perturbations 

In this section we apply the quantum theory of cosmological perturbations on the quan- 
tum, pre-infiationary space-time to extend the study of cosmic inhomogeneities all the way 
back to the Planck era. In addition to the conceptual completion provided by the inclusion 
of Planck scale physics, the resulting framework opens an exciting avenue to extend obser- 
vations into the Planck regime. Before entering into technical details, we summarize here 
the physical idea behind this possibility. 

It is known since the seminal work by Parker in the 60's [85, 86], that a dynamical 
expansion of the universe is able to excite quanta, or 'particles', of test fields out from an 
initially vacuum state. This phenomenon of particle creation is one of the main features of 
QFT in curved space-times, and plays a key role in black hole thermal radiance and in the 
generation of cosmic inhomogeneities during inflation. If k represents a co-moving Fourier 
mode of a test scalar field in FLRW, excitations on that mode may be created if the energy 
scale provided by the space-time scalar curvature is comparable to the physical wavelength 
A = 27ca/k at some time during the evolution. The amount of quanta created during a period 
of expansion in each mode depends on the details of the scale factor a{t) as a function of time. 
Let us focus on the finite range of momenta that is accessible in cosmological observations. 
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The previous argument tell us that, even if those modes are 'born' in the ground state at 
time of the bounce, particles may be created during the evolution. The resulting state, 
e.g. at the onset of inflation, would then depart from the vacuum state at that time as a 
consequence of the non-trivial evolution, and the spectrum of particles created will carry 
information about the pre-inflationary space-time geometry. 

Furthermore, it has been shown in the context of inflation that the predictions for the 
CMB and the distribution of galaxies are sensitive to the details of the state describing per- 
turbations at the onset of inflation [88-91], and concrete observation have been proposed that 
could reveal information about that state [92-94]. In other words, those observations may 
reveal information about the propagation of perturbations before inflation, when quantum 
gravity corrections dominate. 

In the inflationary scenario observable modes have wavelength much smaller than the 
radius of curvature at the onset of inflation (in the cosmological argot, modes are deeply 
inside the Hubble radius). The sometimes implicit assumption in inflationary physics is 
that, whatever happened before inflation, wavelength of interest were much smaller than 
the radius of curvature at any time before inflation. Under this assumption, pre-inflationary 
dynamics for those modes is indistinguishable from an evolution in Minkowski space-time, 
and the use of a vacuum state is justifled. The relevant question is then: is this assumption 
accurate in the pre-inflationary background provided by LQC? More explicitly, consider 
modes with physical wavelength smaller that the radius of curvature at the beginning of 
inflation, and propagate them backward in time until the bounce. Do those wavelength 
generically remain smaller that the radius of curvature of the dressed metric gab during the 
entire pre-inflationary evolution? The detailed analysis of [76, 78] shows that the answer to 
this question is in the negative (sec Fig. 2). While short enough wavelengths (large enough 
momenta) remain always smaller that the curvature radius, there are modes which at some 
time during the evolution have physical size comparable to it. The evolution of those modes 
is sensitive to the space-time curvature and the quantum state at the onset of inflation will 
depart from the vacuum. 

Notice that in LQC the maximum value of the curvature takes place at the bounce 
time and this value is universal, fixed by the quantum geometry and independent of the 
form of the scalar field potential. If we call ku the co-moving scale associated with this 
maximum value of the curvature, we expect excitations with k ^ kp> to be created during 
the evolution, concretely in the Planck regime near the bounce. On the other hand, for 
modes with k ^ k^ pre-infiationary dynamics has negligible effect. From this qualitative 
discussion we may expect observable effects from Planck scale physics in CMB and large 
scale structure if observations are accessible to modes k around or smaller than the universal 
scale kfi provided by LQC. 

In the remainder of this section we provide precise computations that support this quali- 
tative physical picture. We start by specifying the initial condition for both background and 
perturbations at the bounce. We then evolve those perturbations until the end of slow-roll 
infiation, compute the resulting quantum state and the power spectrum for scalar and tensor 
perturbations, and study under what set of initial conditions quantum gravity corrections 
may be sizeable for observable modes. 



41 




FIG. 2: This plot shows: i) The scalar curvature of the effective geometry (red solid line), ii) The 
physical momentum squared (A;/a(t))^, for k = 6 (dotted black line), and A; = 10 (dashed black 
line), and iii) {kji/d{t))'^ , where kji is the co- moving scale associated with the maximum value of 
the curvature (dotted-dashed green line); as a function of cosmic time t. By convention, we choose 
the scale factor of the effective geometry to be one at the bounce, a(0) = 1. Both axes are in 
Planck units. Curvature attains the maximum value at the bounce and decreases very fast after it. 
Modes with momentum k larger than the scale of curvature at the bounce, k > kpt, have physical 
momentum larger than the curvature during the entire evolution (dashed black line). Those modes 
do not 'feel' the curvature and evolve as if they were in Minkowski space-time. On the other hand, 
modes that at the bounce have physical momentum smaller that the curvature, k < kji, quickly 
evolve to become of the same order as the curvature scale (black dotted line), and therefore their 
evolution will differ considerably from that in flat space. At later times those modes also become 
two energetic to feel the space-time curvature. 



1. Initial Conditions 

In the standard inflationary paradigm one specifies 'initial data' for the background and 
perturbations at the onset of slow-roll. From a fundamental point of view, it would be more 
satisfactory to impose initial conditions at the 'beginning' rather than at an intermediate 
time in the evolution of the universe. In classical cosmology the 'beginning' is the big bang 
singularity, and it is not possible to unambiguously defined initial condition at that time. In 
LQC the big bang is replaced by a quantum bounce where physical quantities do not blow 
up, providing a preferred time to specify initial data. 

In the test field approximation, the total wave function naturally decomposes as a 
product \1/ = \&o ® '0; and this form holds as long as back-reaction of perturbations remains 
negligible. We need therefore to specify initial data for both, \[^o and ip. 

• Background. For computational purposes, it is convenient to make the following further 
simplification on the background dynamics. As described in section II A, the background 
wave function \&o can be chosen to be highly peaked along the entire evolution, including 
the deep Planck regime. The 'peak' of that wave function describes an effective geometry 
characterized by the scale factor a(0) = (a(0)), which satisfies the effective equation (2.58). 
Because the dispersion of \l/o remains very small during evolution, it is convenient to ignore 
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quantum fluctuations in our computations, by making a 'mean field' approximation in which 
the expectation values of powers of background operators, such as a and Ho^ are replaced 
by the same powers of their expectation. For instance, in the evolution of quantum inho- 
mogeneities given by Eq. (3.11), this is equivalent to replace a ~ a. At the practical level 
this is an excellent approximation, e.g. numerical errors in simulations turn out to be larger 
than those introduced by the mean field approximation. 

In subsection IV A we described the effective pre-inflationary background arising in LQC 
for the representative example of a quadratic potential. In that effective geometry initial 
data is entirely specified by the value of the scalar field at the bounce, 0b, and, unless 0b 
lies in a small region R around 0b=O, the evolution generically finds an infiationary phase 
at late times compatible with WMAP observations [77]. Therefore, we will choose ^E'o to be 
a state sharply peaked in an effective trajectory specified by a value oi 4>b that lies outside 
the region R. 

The effect of choosing different values of 0b can be understood using the effective 
equations (2.58) together with numerical simulations. On the one hand, immediately after 
the bounce the background evolution is entirely dominated by quantum gravity effects, 
and it is largely insensitive to the concrete value of 0b- Except for very small momenta 
k, the times at which perturbations Qk and Tk 'feel' the space-time curvature is precisely 
just after the bounce (see Fig. 2). Therefore, the features that those modes acquire during 
the evolution turn out to be quite insensitive to the value of 0b- On the other hand, 
different values of 0b do modify significantly the space-time geometry at later times. 
The larger 0b, the longer it takes to reach the end of slow-roll inflation, or, equivalently, 
the larger the amount of expansion of the universe between the bounce and the end of 
slow-roll. A larger amount of expansion implies that observable modes had larger physical 
momentum at the time of the bounce. Because by convention we fix the scale factor at 
the bounce ob = 1 (rather than atoday = 1), the effect of choosing different values of 0b 
essentially translates into a change in the range of co-moving momenta k relevant for 
observations, moving to larger fc's as 0b increases. If [/Cmmj ^max ~ 2000/Cmin] is the window 
covered by WMAP, we have, for instance, A';j„in ~ 2.8 x 10^"^ for 0b = 1, ^ 0.14 for 
0B = 1-1 and /Cmin ~ 8.2 for 0b = 1.2. The physical momentum /c/otoday of modes observed 
today is of course the same in all cases, but the convention qb — 1 makes that different 
amount of expansion (i.e. different 0b) translates into different co-moving k for those modes. 

• Perturbations. As already occurs in classical space-times, quantum fields in quantum 
cosmological backgrounds does not admit a preferred state that we can call the vacuum. In 
backgrounds with large enough number of isometrics, e.g. Minkowski or de Sitter space- 
time, a preferred ground state can be singled out by imposing symmetry in combination 
with regularity conditions. In our quantum FLRW we follow the same criteria, and look 
for quantum states ip invariant under the isometrics of the background, spatial translations 
and rotations, with appropriate ultraviolet behavior. In section III B 2 we summarized the 
construction of the Hilbert space T-ti of 4th-order adiabatic states. In "Hi, the family of 
4th-order adiabatic vacua is the preferred set of initial conditions selected by symmetry 
and regularity requirements. This is the set of initial data we choose for perturbations. 
As opposed to Poincare or de Sitter invariance, symmetry under spatial translations and 
rotations is not restrictive enough to select a unique state, but it substantially narrows 
down the possibilities. The next subsection will summarize the time evolution of different 
choices of initial state within the family of 4th-order adiabatic vacua, and will show that 
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quantities of interest such as the power spectrum of observable modes, are all very similar. 

Physically, the choice of a 4th-order adiabatic vacuum at the time of the bounce corre- 
sponds to assume 'initial quantum homogeneity'. One is requiring that the portion of the 
universe corresponding to our observable patch at the time of the bounce is as homogeneous 
as quantum mechanics allows, i.e. only vacuum fluctuation of inhomogeneities are present. 
This is a strong assumption. The motivation comes from [76, 78]: 

• In a universe containing a phase of inflation lasting at least for 60 e-folds, the physical 
size of observable universe was very small at the bounce time, < 10£pi, for the solutions 
of interest. 

• The 'quantum degeneracy force' responsible of the bounce has a diluting effect that 
may produce homogeneity at scales of the order of the Planck length at the bounce. 
This is the new ingredient that LQC provides at the time of the bounce to produce 
homogeneity at Planck scale distances. 

• There is a precise sense in which the assumption of quantum homogeneity captures a 
quantum version of the Weyl curvature hypothesis [95]. 

2. Power Spectrum 

Our task is to use the equations of the quantum theory summarized in section III B to 
compute the state of cosmic inhomogeneities at the end of the inflationary epoch, by starting 
from the initial condition specified above for background and perturbations at the time of 
the bounce. 

Due to computational limitations, it is convenient to restrict numerical simulations to 

backgrounds for which the bounce is kinetic energy dominated, where it has been shown 
that quantum fluctuations of remain very small along the entire evolution. Several 
numerical simulations have been carried out for effective backgrounds with initial conditions 
05 e (0.93, 1.5), which turns out to be the most interesting range [76]. It is not expected 
that new features appear for larger values of 0b, but computational limitations make difficult 
to check it explicitly. 

For perturbations, simulations have been carried out using different choices of 4th-order 
adiabatic vacua, and the results are all very similar. Fig. 3 and 4 are obtained by using 
the 'obvious' or 'standard' 4th-order vacuum at the bounce time fjB (see [4] for precise 

definition), and they show the relevant information of the evolved state. 

First of all, to gain intuition on the effect of the pre-inflationary evolution, we compare 
the evolved state with the natural vacuum during inflation, the so-called Bunch-Davies 

— * 

(BD) vacuum. Fig. 3 shows the number Uk of 'excitations/particles' with momentum k per 
comoving unit volume in space and momentum, contained in the evolved state relative to 
the BD vacuum during inflation. The plot is computed for 0^ = 1-15 but, as explained in 
subsection IV B 1, it is not altered by choosing a different value inside our family. Changing 
the value of 0b has essentially the effect of shifting the location of the observationally relevant 
window [/cmin, ^max ~ SOOO/cmin] in the horizoutal axes of the plot, which moves steadily to the 
right as 05 increases. Fig. 3 is in good agreement with the qualitative arguments presented 
at the beginning of section IV B. Namely, the pre-inflationary evolution affects modes with 
low k, for which a considerable amount of excitations have been 'created'. On the contrary, 
modes with large k remain in the ground state at the onset of inflation. As it was expected. 
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FIG. 3: Number Uk of scalar 'excitations/particles' with comoving momentum k in the interval 
[k, k + dk], per comoving unit volume contained in the evolved state as compared to the BD vacuum 
during inflation. The plot is computed for (ps = 1-15 and for the 'obvious' 4th-order adiabatic 
vacuum at the bounce. The horizontal axes is in Planck units. 
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FIG. 4: Ratio of the LQC power spectrum for scalar perturbation to the standard inflationary 
power spectrum. Crosses show the ratio for different values of k. The LQC power spectrum 
oscillates rapidly for small k. The solid curve averages over bins of width AA; = 0.5. The inset 
shows a zoom-in of the interesting region around k=9. 



for k > ^ 7.7 (recall that kji is the comoving scale associated with the scalar curvature 
of the effective metric at the bounce), the number of BD particles contained in the evolved 
state is very close to zero. Therefore, if fcmin ^ kn, that corresponds to 0b > 1.2, the evolved 
state is indistinguishable from the BD vacuum for observable modes. For (pB ^ 1-2 the state 
at the onset of inflation differs significantly from the vacuum for modes in the interesting 
window and, as analyzed in detail in [88-91], those deviations have an important effect on 
the predictions of inflation for the spectrum of cosmic inhomogeneities, specially regarding 
non-Gaussianity. There exist concrete proposals for observables in the CMB [93, 94] and 
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in the distribution of galaxies [92, 93] that should be sensitive to the effects of the created 
particles. 

A quantity of direct observational interest is the power spectrum of tensor and scalar 
perturbations 



where all quantities are evaluated at the end of inflation, H is the Hubble rate, and qk{t) and 
ek{t) are the Fourier modes of scalar and tensor perturbations, respectively. Fig. 4 shows the 
relation between the LQC power spectrum computed with the evolved state and the standard 
inflationary power spectrum that assumes the BD-vacuum, for scalar perturbations. The 
conclusions are similar to the ones obtained from Fig. 3, namely for 0^ ^ 1-2 the power 
spectrum of observable modes is indistinguishable from the standard inflationary predictions. 
For smaller values of 0b deviations become sizable for modes of observational interest. For 
instance, for (})b = 1-15 we have fcmm ~ 1 and deviations from standard prediction will 
appear for modes with £ < 30 in the WMAP angular decomposition. These deviations are 
inside current uncertainties. However, the fact that the state for perturbations differs from 
the BD-vacuum opens a window to observe those effects. 

The analogous plot for tensor modes has the same form as Fig. 4, and the conclusions are 
also the same [76] . In particular, there are no important corrections for the tcnsor-to-scalar 
ratio, although the inflationary consistency relation, which relates the tensor-to-scalar-ratio 
and the tensor spectral index, is modified [76] . 



The last step is to check whether there exist a big enough set of physical states ip on 
the Hilbert space for which the truncation underlying our quantum theory, the test field 
approximation, holds during the entire evolution. This is an intricate question because: 
i) It requires a detailed analj^ical control of the necessary regularization on states and 
composite operators on our Hilbert space; ii) Numerical implementation of those techniques 
are necessary to check self-consistency at any time during the evolution^ dealing with the 
subtleties of having numerical control on the subtraction of quantities that tend rapidly to 
infinity, during a period that covers around 11 orders of magnitude in energy density. 

Section III B summarized the necessary tools to check self-consistency and pointed out 
that a sufficient condition is that the energy density in perturbations (p) be negligible 
compared to the background (po) o-t (iny time during the evolution. Fig. 5 shows the result of 
the numerical evolution of the energy density for scalar perturbations (analogous results hold 
for tensor perturbations). The plot shows the ratio (pq) / (po) for a background corresponding 
to (f) — 1.23 and the 'obvious' 4th adiabatic order vacuum specified at the bounce. This ratio 
remains small for the entire evolution, including the Planck regime. The initial condition 
= 1.23 corresponds to /cmin ~ 30, therefore the number of excitations over the BD state 
on observable modes is negligible (see Fig. 3) for this background. Additionally, there exist 
an analytical argument [76] ensuring that, given a state for perturbations for which back- 
reaction is negligible, there exist a well defined neighborhood of that state with the same 
property. Each of those provide a state at the beginning of slow-roll indistinguishable from 
the BD vacuum. They provide therefore, viable extensions of the standard infiationary 
scenario that includes Planck scale physics [76, 78]. 



^3 

Pr{k) = h—^\ek 




(4.3) 



3. Self- consistency 
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FIG. 5: Ratio of the energy density of scalar perturbation to the background energy density as a 
function of cosmic time. The initial conditions were chosen as (pB = 1-23 for the background, and 
the 'obvious' 4th-order adiabatic vacuum at the bounce for perturbations. Slow-roll inflation starts 
about 3 X 10^ Planck seconds after the bounce. During the entire evolution the ratio remains small. 
This example constitutes a self-consistent extension of the evolution of cosmic inhomogeneities to 
include the Planck era. 

For the range (pB < 1-2 there are only upper bounds for {pg) which are far from being 
optimal. At the present time there are no explicit computations for which the test field 
approximation is satisfied for (pB in that window, and additional work is required to establish 
the self-consistency of our truncation scheme. 

V. CONCLUSIONS 

One of the most pressing questions a quantum theory of gravity has pertains to both 
theoretical and observational issues in cosmology. In the theoretical front the standard 
model is based on General Relativity that possesses an initial singularity, a signal that the 
theory breaks down at some point. On the observational front, the CMB spectrum poses 
very stringent conditions for any theory of the early universe. One of such scenarios is given 
by the inflationary paradigm, that explains very successfully the detailed structure of the 
inhomogeneities seen in the CMB as an imprint of quantum fluctuations of certain fields 
just before the inflationary phase. Can one have a formalism that provides a satisfactory, 
nonsingular description both at the Planck scale and at the onset of inflation? Interestingly, 
loop quantum cosmology allows to answer both questions in the affirmative. 

As we have described in this Chapter, when one considers the homogeneous degrees 
of freedom, the so called 'background geometry', the formalism provides precise singularity 
resolution, replacing the classical big bang with a big bounce. The dynamics of semiclassical 
states is very well described by an effective theory that captures the leading quantum gravity 
effects and allows one to describe the spacetime geometry in terms of an effective background 
metric. 

The inflationary scenario is very powerful to explain in great detail many features of the 
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observed CMB spectrum. It is however, incomplete in various directions. In particular, it is 
based on General Relativity where the spacetimes under consideration are past incomplete, 
that is, singular. As we have described in detail, one can indeed extend the scenario back 
in time to the Planck scale. For that one needs two new ingredients. The first one is a for- 
malism that allows one to treat quantum perturbations of the spacetime metric propagating 
not on a classical spacetime, but rather on a quantum spacetime. The second ingredient 
involves consistency conditions that ensure us that one can 'evolve' the quantum perturba- 
tions back to the Planck scale without violating the approximations that yield vahdity to 
the formalism. As we have seen one can indeed consistently consider the extension of the 
inflationary scenario. 

Perhaps the most important question is whether this extension to the quantum bounce 
provides a window for Planck scale physics to be observed in the CMB. As we have described, 
the sector of the parameter space that has been explored provides predictions that are fully 
consistent with the standard inflationary scenario, under current observations. Further 
explorations are needed to decide whether the scenario provided by LQC is both consistent 
in the full parameter space and provides us with distinct testable predictions. 
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